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This paper reviews recent advances in continuous-time quantum walks (CTQW) and their application to 
transport in various systems. The introduction gives a brief survey of the historical background of CTQW. After 
a short outline of the theoretical ideas behind CTQW and of its relation to classical continuous-time random 
walks (CTRW) in Sec. 2, implications for the efficiency of the transport are presented in Sec. 3. The fourth 
section gives an overview of different types of networks on which CTQW have been studied so far Extensions 
of CTQW to systems with long-range interactions and with static disorder are discussed in section V. Systems 
with traps, i.e., systems in which the walker's probability to remain inside the system is not conserved, are 
presented in section IV. Relations to similar approaches to the transport are studied in section VII. The paper 
closes with an outlook on possible future directions. 
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I. INTRODUCTION 

Transport of mass, charge or energy is the basis of many 
physical, chemical or biological processes. Such transfer 
mechanisms and their efficiency depend on the underlying 
structure of the system, range from polymer physics to solid 
state physics to biological physics to even quantum compu- 
tation, see Refs. |l-6] for reviews. The underlying structures 
could be, for example, simple crystals, as in solid state physics 
I?!], more complex molecular aggregates like polymers [SJ, or 
general network structures 

In quantum mechanics, the potential a particle is moving 
in specifies the Hamiltonian of the system, which determines 
the time evolution. For instance, the dynamics of electrons in 
a simple crystal is described by the Bloch ansatz iItI l9l [Toll . 



2 



which mirrors their behavior in metals quite accurately. Over 
the years these models have been refined and extended to ad- 
dress different phenomena, such as the dynamics of atoms 
in optical lattices and the Anderson localization in systems 
with energetic disorder [11.1 . In quantum chemistry, Hiickel's 
molecular-orbital theory allows to define a Hamiltonian for 
more complex structures, such as molecules fl?]. This is also 
related to transport in polymers, where the connectivity of the 
polymer plays a fundamental role in its dynamical and relax- 
ational properties 1 13]. There, (classical) transport processes 
can be modeled by continuous-time random walk (CTRW) ap- 
proaches liSt j6l [st] . An important example is the motion of 
Frenkel excitons, whose high-temperature dynamics is often 
governed by a master equation with an appropriate (classical) 
transfer operator which determines the temporal evolution of 
the system ifTilfTsll . CTRW have proven to be a very useful 
and general tool in describing incoherent transport in various 
settings, ranging, for instance, from percolation on fractals 
fl6] to electronic energy transfer in glassy systems |17, 1^. 
Also, anomalous diffusional behavior has been successfully 
modeled by CTRW using suitable waiting time distributions 

MM- 



A. Quantum Walks 

The tight-binding approximation used in sold state physics 
as well as in Hiickel's theory is equivalent to the so-called 
quantum walks, which model purely coherent quantum dy- 
namics of excitations on networks P, '2l'-'27l. What mat- 
ters is that the constituting elements (spins, atoms, molecules, 
etc.) are of the same type, in the simplest cases acting as two- 
level systems. Now, the extension of classical random walks 
ideas to the quantum domain is not unique and allows vari- 
ants, of which two types capture most of the interest: discrete- 
time quantum (random) walks (introduced by Aharonov et ai, 
using an additional internal "coin" degree of freedom ll22ll ). 
and continuous-time quantum walks (CTQW) (introduced by 
Farhi and Gutmann, where the connection to CTRW uses 
the analogy between the quantum mechanical Hamiltonian 
and the classical transfer matrix [2Jj]). Recently, Strauch has 
shown how these two versions are related Q. Thus, CTQW 
and CTRW act as the two extreme cases of purely coherent 
and purely incoherent transport, respectively. 

In quantum information theory, quantum walks are used ex- 
tensively as algorithmic tools for quantum computation ^m■ 
Here, the most prominent examples are Shor's algorithm 127|] 
and Grover's algorithm [28]. The latter is a quantum algo- 
rithm for finding an item in an (unsorted) database of qubits. 
Such an algorithm can also be related to CTQW |29, 30]. The 
idea of using computers built on quantum mechanical prin- 
ciples instead of classical computers goes back to Feynman, 
who already formulated an early version of quantum walks in 
1982 I3lll . From a conceptual point of view, quantum walks 
are closely related to quantum cellular automata, see, e.g., 
ll32ll . Using CTQW, the transfer of information in complex 
systems was discussed by Christandl et al. in the context of 
perfect state transfer [331 ■ 



CTQW are also in close relation to the so-called quantum 
graphs (QG), see, e.g., fPUSTtj. However, unlike CTQW, 
QG take into account detailed properties of each bond of the 
graph explicitly; in QG, bonds may be directed and have dis- 
tinct lengths. The connections between discrete Laplacians on 
discrete QG and periodic orbits were recently investigated by 
Smilansky S. 



B. Experimental implementations of CTQW 

Simple theoretical models have always been very useful 
for our understanding of physics. In quantum mechanics, 
next to the harmonic oscillator, the particle in a box provides 
much insight into the quantum world (e.g. l39tl ). The problem 
of a quantum mechanical particle moving in an infinite box 
has been reexamined in [40-42]: Remarkably, even this sim- 
ple system shows complex but regular spacetime probability 
structures (so-called quantum carpets). In solid state physics 
and quantum information theory, one of the most simple sys- 
tems is associated with a particle moving in a regular periodic 
potential. 

In recent years, arrangements close to ideal theoretical 
models have been tailored using, e.g., ultra-cold atoms in opti- 
cal lattices, see ||43|] and references therein. For both quantum 
walk variants, experimental implementations have been pro- 
posed, based on microwave cavities 144.1 . on Rydberg atoms 
i45ll or on ground state atoms ||46|] in optical lattices or in 
optical cavities l47ll . or using the orbital angular momentum 
of photons ll48ll . Other experimental proposals connected to 
CTQW are based on waveguide arrays ll49ll or on structured 
clouds of Rydberg atoms IdOII . 



C. Transport in complex systems 

The range of systems where CTQW can be used to study 
transport is clearly not restricted to simple (symmetric) mod- 
els. In particular, CTQW have been proven to be very useful 
in describing the dynamics of excitations in various complex 
systems, as will be shown throughout this review. 

Since CTQW can be related to the tight-binding approxima- 
tion in solid state physics, CTQW are also related to chains of 
coupled spin 1/2 particles. If all but one spins are prepared in 
the same state, say, spin down, and a single spin in the spin 
up state, one can map this onto a chain of coupled two-level 
systems with a single initial excitation at the node which is 
identified with the spin up particle ll24ll5]l l52ll . The essential 
dynamics is then equivalent to a CTQW on a linear network. 

However, also the dynamics in topological disordered sys- 
tems can be modeled by CTQW. Take for instance a gas of 
highly excited ultra-cold Rydberg atoms. At ultra-low tem- 
peratures the configuration of the atoms can be thought of as 
being frozen on the time-scale on which the following trans- 
port process is taking place ifssl - lssll : The dynamics is started 
by exciting one (or several) of the Rydberg atoms in a higher 
state which is resonantly coupled to a lower Rydberg state 
ll45il50}l . Since essentially only two different Rydberg states 
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are involved, one can map this again onto a network of cou- 
pled two-level systems. 

Also biological systems share some properties of simple 
network models. Recent experiments on the light-harvesting 
complexes of algae have renewed the interest in studying the 
dynamics of excitations in such systems ll56i IstIi . Here, an 
initial excitation (a Frenkel exciton) is created by absorbing 
a solar photon. The exciton is then transfered along a net- 
work of (bacterio-) chlorophylls (BChl), the chromophores, 
to the reaction center, where the excitation energy is con- 
verted into chemical energy. Now, the network of the BChls 
can be considered to be static and stable with defined cou- 
plings between the BChls, at all (even ambient) temperatures 
of interest. Moreover, the Frenkel exciton can be viewed as 
a quasi particle moving along the network. Therefore, also 
here the transport can be modeled by CTQW, at least at (very) 
low temperatures. The incoherent exciton transport in den- 
drimers can be efficiently modelled by random walks, see, 
for instance, |58-61]. However, for certain systems there 
is also experimental evidence for coherent interchromophore 
transport processes ll62i 163^1. Recent investigations of biologi- 
cal light-harvesting systems range from the Fenna-Matthews- 
Olsen complex p6] to marine algae jSTll . There have also 
been recent theoretical efforts to understand the coherent fea- 
tures of exciton dynamics found in recent experiments, see, 
e.g., iBMsfl. 

Clearly, such biological systems cannot be thought of as be- 
ing isolated from their environment. Therefore, in these cases 
the purely coherent description by CTQW is of limited value. 
However, it is also possible to extend the CTQW approach to 
take the coupling to the environment into account, which leads 
to models related to quantum master equations Ii69l1 . see also 
Sec lylTB] 



stance, is given by 
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Take as an example a ring-like network of N nodes, with 
only nearest neighbor connections and + 1) = then 
the matrix A reads: 



/2 -1 ... -1\ 
' -1 2 -1 ... 







V-i 



-1 2 -1 
-1 2/ 



(3) 



Writing A in a quantum mechanical fashion using the projec- 
tion operators \l) {k\, leads to 



A = J22\1}{1\-\1-1}{1\-\1 + 1){1\ 



(4) 



Having now specified the network through the connectivity 
matrix A, one is interested in the dynamics over the network. 
First, one distinguishes between purely classical, i.e., incoher- 
ent dynamics and purely quantum-mechanical, i.e., coherent 
dynamics. The main question one wishes to answer is: What 
is the probability to be at node k after time t when starting at 
node j? While here the focus is mainly on localized initial 
conditions at a single node j, in general, the initial state can 
be any distribution involving all the nodes. 



II. DEFINITIONS AND TERMINOLOGY OF CTQW 
A. Connectivity of networii 

Networks involved in the CTQW and CTRW dynamics are 
characterized by the form in which their sites are connected. 
Starting point of a network is a collection of N nodes, which 
are then joined by bonds, the connectivity matrix A mirroring 
the way in which the bonding occurs. 

The N X N connectivity matrix A has the elements Akj, 
where 

{/, for k^j 
— 1 if fc and j connected Akj — Ajk (1) 
else, 

where fj is the number of bonds emanating from node j. This 
matrix has interesting and useful properties: (a) A is real and 
symmetric, (b) all its eigenvalues A„ are real and A„ > 0, and 
(c) A has a single smallest eigenvalue which is Ai — 0. 

One can associate with every node of the network a ba- 
sis vector in an A^-dimensional vector space. Now, these ba- 
sis vectors form a complete orthonormal basis, which, for in- 



B. Transfer matrix for CTRW 

We start by considering classical transport. Let the initial 
node be j, such that the initial state of the system is \j), and 
denote the transition probability to go in time t from node j to 
node k by pkj{t). Therefore, the initial condition is {k\j) = 
Pfc.j(O) = ^k.j, where Si^ j — 1 for k = j and otherwise. 
The state after time t is \ j; t), such that the overlap with node 
fc reads {k\j;t) =pk^j{t). 

The dynamics resulting in the state \ j; t) follows from the 
transition rates per unit time between two nodes. Those tran- 
sition rates are the elements of the so-called transfer matrix 
T, Tfcj = (fc|r|j), which is therefore related to the spatial 
gain or loss. One can now write the temporal change of the 
probability after time e <C 1 as 

Pk,] (e) = Pfe J (0) + eTfcj = 4 J + eTkj . (5) 

Assuming a Markovian process, the following master equa- 
tion can be shown to hold [5, 6]: 

Y,P'^At)^J2^kiPiAt)- (6) 
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This equation defines the CTRW, which (depending on the 
specific form of the transfer matrix) have been appHed to var- 
ious problems in physics, chemistry, biology, and also social 
sciences. 

In the simplest case, where the rates for all bonds are equal, 
say, 7, the transfer matrix is related to the connectivity matrix 
through T = —jA. Then the master equation describes a 
diffusive motion over the network. 

The formal solution of Eq.© is Pk,j{t) = (/cje-^^'lj) = 
(fc|e~^^*|j'). Denoting the eigenstates of A by |q„), one has 



Pk.j{t) = {k\e 



-7At| -N _ V^, 



\J) = 



\qn){qn\j). (7) 



Since the eigenvalues are positive (A„ > for 7i > 1 and 
Ai ~ 0), the long-time limit follows directly: For t ^ 1 
in the sum of Eq.(|7]i all exponential terms but one decay 
quickly to zero. The only term which survives is the one for 
Ai = 0, with the corresponding eigenstate — j^J^i 10- 
Therefore, the long-time limit of all transition probabilities is 
limt^oo Pkj {t) — l/N, which is independent of the connec- 
tivity of the network. This means that every CTRW whose 
transfer matrix follows directly from the connectivity matrix 
will eventually decay at long times to the equipartition value 
a fact which is sometimes referred to as "ground state 
dominance". 



C. Hamiltonian for CTQW 



extremes of transport on the same topology. On the one hand 
there is the purely incoherent CTRW, while on the other hand, 
there is Schrodinger's equation which now defines CTQW. 
Assuming again that all transition rates between different con- 
nected nodes are the same, the Hamiltonian can be directly 
related to the connectivity, H = jA. Therefore, the underly- 
ing topological network is determined both for CTRW and for 
CTQW by the connectivity matrix A. This allows to study the 
role of the connectivity in parallel for CTRW and for CTQW. 

Let us donote the eigenvalues of H by £'„ and the eigen- 
states of H by l^'n). Evidently, the eigenvalues and eigen- 
states of H and T are practically the same; nonetheless, we 
will keep the distinction between H and T due to later pur- 
poses. Now, the quantum-mechanical transition probabilities 
read 



(10) 



Unlike the situation for CTRW, for CTQW there is no unique 
long-time limit of TTkj (t) due to the unitary time evolution. In 
order to compare to the long-time behavior of the CTRW one 
uses the long time average 
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Xk,j = lim - / dtiTk.j{t) 



T 



|*n)(*„|i)(j|*m)(*m|fc), (11) 



Turning now to the quantum-mechanical dynamics, one 
faces the fact that the transport has to be formulated in Hilbert 
space. For this one assumes that the states \i) representing the 
nodes span the whole accessible Hilbert space. As before, it 
is assumed that the states are orthonormal and complete, i.e., 
= 'Y^- \ — 1. The dynamics is then governed 
by a specific Hamiltonian H, such that Schrodinger's equa- 
tion for the transition amplitudes akj (t) = {k\j; t) reads 



d 



':'^Hk,i ai.j{t). 



(8) 



Similar to the CTRW, the formal solution for the transition 
amplitudes is given by 



-iHt 



(9) 



where e is the quantum-mechanical time-evolution oper- 
ator. The transition probabilities follow as usual as tt^ j (t) = 

Comparing Eq.® to Eq.® one immediatly notices the 
very similar structure of the two equations, except for the 
imaginary unit i appearing in Eq.®. However, while Eq.© 
is an equation for the transition probabilities pkj (t), Eq.® is 
an equation for the transition amplitudes akj (t). 

One can push the similarities further by identifying the 
quantum-mechanical Hamiltonian H with the classical trans- 
fer matrix T, i.e., H = —T. This approach, pioneered by 
Farhi and Gutmann in Ref. ll2lll . allows to compare the two 



where 5E„,Er,^ = 1 if En — E,n and otherwise. The long- 
time average Xkj still depends on the initial and final nodes j 
and k. 



D. Example: A discrete ring 

As an example which illustrates the differences and sim- 
ilarities between CTQW and CTRW we consider a discrete 
ring of N nodes, whose connectivity matrix is given above 
in Eq.lO. This one-dimensional structure with the periodic 
boundary conditions jA^ + 1) — |1) allows for a full analyti- 
cal solution. When realizing that the matrix A is nothing but 
a tight-binding matrix for a particle moving in a regular one- 
dimensional crystal, one readily obtains the eigenvalues and 
eigenstates from a Bloch ansatz |0[T3l- Namely, the Bloch 
states are linear combinations of the localized states | j) and 
are given by 



1 ^ 



(12) 



Now the energy is obtained as 

Eg = 2 - 2cos^ 

For small 6 the energy is given by Eg 
the energy spectrum of a free particle. 



(13) 

9'^ which resembles 
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By inverting Eq.(fT2]i one may describe the state \j) local- 
ized at node j as a Wannier state [i 



(a) N=21 



(a) N=20 



(14) 



Since the states \j) span the whole accessible Hilbert space, 
one has {k\j) = 6kj and, via Eq.(fT2|l. also (5'e'|$6i) — Se'O- 
Then the transition amplitude reads 17211 



e,0' 

l_J2^-^Eet^-^eik^n, (15) 



In an analogous way one obtains the classical transition prob- 
abilities [,7Z1 



-iSik-j) 



(16) 



because the Bloch states are also the eigenstates \qn)- 

The periodic boundary condition for aid lattice of size N 
requires that 6 = 2mr/N with n e]0,N]. Now Eq.([T5ll is 
given by 



ajkit) 



i2t cos(2mT /N) -i2Tvn(k-j)/N 



N ^ 



(17) 



For small Q, this is directly related to the results obtained for 
a quantum particle in a box [40-42], because then £"„ ^ n^. 
Indeed, some features found for the particle in a box can also 
be found in the case of a CTQW on a discrete ring. For the 
particle in a box the initial condition is restored after some 
revival time. In fact the probability to find the particle at a 
certain position in the box is a periodic function. 

Analogously but not completely similarily, for the CTQW 
on the ring, the initial (localized) condition is only partially 
restored. Only small rings of sizes iV = 1,2,3,4, and 6 
lead to a full revival of the initial condition. All other sizes 
only lead to partial revivals. The revival time r is given by 
(^k,j{j) — afc,j"(0). Since the transition amplitudes are given 
as a sum over all modes n, see Eq.dTTll. one has for each mode 
n its revival time Ii72il 



[1 + cot2(ri7r/A^)], 



(18) 



where r S N (without any loss of generality one sets r — 1). 
From Eq.jTSll one finds that r„ > t„_|_i for n e]0, A^/2] and 
Tn < TnJ^x for ji g] A^/2, N\. For certain values of n, Tn will 
be of order unity, e.g. for n — N/2 t„ = 7r/2. However, for 
n « N, Eq.([T8]l yields t„ = N^/2nn'^ = Ta/n^, which is 
analogous to the particle in the box and where tq is a universal 
revival time. 

Because the revival times t„ have large variations in value, 
one compares these times to the actual time needed by the 
CTQW for travelling through the lattice. As mentioned ear- 
lier, interference effects in the return probability 7ri.i(i) are 




node k 



node k 



FIG. 1: Contour plot of the probability for a CTQW on a circle of 
length (a) A'^ = 21 and (b) = 20 over long times t. Dark regions 
denote high probabilities. From i72h . 



seen after a time t w N/2. The first revival time has to 
be larger than this, because there cannot be any revival un- 
less the wave reaches its starting node again. Our calcula- 
tions ll72ll suggest that the first revival time will be of order 
tq. From Fig.[T]one sees that the first (incomplete) revival oc- 
curs for iV = 20 at t « 70 > 20V27r and for = 21 at 
t « 75 > 2lV27r. 

Since the transition probabilities are easily obtained from 
the Bloch ansatz, also the long-time averages can be estimated 
analytically. Depending on the number of nodes in the net- 
work, the long-time averages are, [|72ll . for even N: 



Xk,j 



and for odd N: 



1/N^ k^j 

1/2N k = j andk = j + N/2 



Xk,j 



(19) 




(20) 



Thus, unlike in the classical CTRW case the long-time average 
is not equipartitioned but "remembers" the initial condition. 
Moreover, also the symmetry arising from the even or odd 
number N of nodes is reflected in the Xk,j ■ 



III. EFFICIENCY OF CTQW AND CTRW 

The performance of CTQW and CTRW depends to a large 
extent on the connectivity of the underlying network, i.e., on 
its topology. One focuses on the return probabilities ttj^ (t) 
and pjj{t). Now, if these probabilities decay quickly with 
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time, the probabilities 1 



(t) and 1 — pj,j (t), to be at any 



but the initial node, increase quickly. This implies a fast trans- 
port through the network. In order to make a global statement 
on the performance, one considers the average return proba- 
bihties 1,73,1 



for CTQW 



and 



(21) 



(22) 



Again, a quick decay of Tf{t) [p{t)] implies a fast propaga- 
tion through the network, while a slow decay implies a slow 
propagation. 

For CTRW the average return probabilities simplify consid- 
erably; one has namely 



n j 

^E' 



-A„t 



(23) 



which only depends on the eigenvalues A„ of the transfer 
matrix T but not on its eigenstates. By making use of the 
Cauchy-Schwarz inequality one gets a similar expression for 
CTQW. The average return probabily 7f (t) is related to the av- 
erage return amplitude a{t) by Il73ll 



(24) 



Similar to p{t) for CTRW, only depends on the eigen- 

values En of the Hamiltonian H: 

= |^E^"'''"*(*"iEi^')oi*'^ 

n j 

- I-E 

\n ^ 



(25) 



Having a quantity which only requires the calculation of the 
eigenvalues considerably shortens the computation time, es- 
pecially for very large systems. 

Clearly, due to the oscillating terms in Eq.(l25]l. the lower 
bound will oscillate in most cases. Thus, to compare 

to the decay of the classical p{t) one uses the envelope of 
|a(t)p. Note that when one identifies the Hamiltonian of the 
CTQW with the transfer matrix of the CTRW, the eigenvalues 
of both are the same, i.e.. En = A„. Thus, formally the differ- 
ence in dynamics is only due to the different functional form 
of Eq.(|23l) and Eq.(|25]). Nonetheless, this can lead to drastic 
effects. 



As already discussed, for CTRW the long-time limit of the 
transition probabilies reaches the equipartition value 1 /N. In 
the same way, the long-time limit of p{t) is given by 1 /N. In 
contrast, for CTQW neither Tf{t) nor decay to a given 

value at long times, but rather oscillate around the correspon- 
nding long-time average which for 7f(i) is given by ll74ll 



1 r 
= lim - / dt7t{t) 



(26) 



Again, on can obtain a lower bound which does not depend on 
the eigenvectors 1 74] : 



X 



> 



Xib- 



(27) 



A. Examples 

The difference in the global efficiencies of CTQW and 
CTRW is illustrated by considering two distinct examples: the 
discrete ring of N nodes and the star with one core node and 
iV — 1 nodes attached to it, see Fig.|2] 

The eigenvalues (and also eigenstates) of the discrete ring 
have already been discussed. From these results it is straight- 
forward to calculate the average return probabilities. Since the 
ring's eigenstates are Bloch states, one can easily verify that 
the lower bound is exact, i.e., in this case also 7f(t) 

does not depend on the eigenstates. Therefore, the average 
return probabilities are given by ItsIi 



I 1 ^ 

W{t) = I X! exp[-i2t(l - cos(27m/A^)) 



n=l 



(28) 



for CTQW and 



1 ^ 

= ]^ E exp[-2t(l - cos(2^n/A^))] (29) 



for CTRW. In the limit of large N, i.e., for continuous 9 = 
2TTn/N, one can replace the sums by integrals which leads to 



and to 



dE p{E) e^p{-Et), 



dE p{E) eyiY>{-iEt) 



(30) 



(31) 



where p{E) is the density of states (DOS). Then one obtains 
for CTQW 7f(t) - |Jo(2t)P, which for i > 1 can be ap- 
proximated by 7f(t) ^ sin^(2i + 7r/4)/f [75]. Since the max- 
imum of the sine function is 1, the envelope of 7f(t) decays 
as . For CTRW there is only a single sum and no ad- 
ditional quadrature, such that a similar calculation leads to 
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p{t) ^ t~^/'^. This temporal decay will also be present for 
finite systems. The time-range over which it will be visible 
depends on the value of N . Obviously, for small N the long- 
time behavior will be reached faster than for larger N . 

Thus, we just established that for CTQW on a ring the ex- 
ponent of the decay power-law is twice as large as the one for 
CTRW. This fact holds for a wide class of systems. Namely, 
all networks whose density of states p{E) = jj J2n=i ^{-^ ~ 
En) follows a power-law, p{E) ~ {EE„^ — E'^Y, will show 
this feature. Here, _E,„ is the maximal eigenvalue (one as- 
sumes the minimum eigenvalue to be zero). Since the in- 
terest is in the large t behavior, p{t) will be mainly deter- 
mined by the small eigenvalues, such that for t 1 one has 
p{E) E" .11 is straighforward to show that for CTRW 



p{t) - 



(32) 



This scaling feature at not too short times is well known, see, 
e.g., ll76ll . where 2(1 + v) = ds is sometimes called the spec- 
tral or fracton dimension. For CTQW, also will be in 
general determined by the small eigenvalues. For p{E) ^ E^ 
one obtains — p{t). Here, all quantum mechanical os- 
cillations vanish, because one considers only the leading term 
of the DOS. Therefore, the envelope of the lower bound for 
CTQW scales as iItI 

^-2(1+1.)^ 



env[|a(i)|" 



(33) 



In other cases of interest (e.g., stars, see below) the DOS has 
highly degenerate eigenvalues |74]. As an extreme case there 
may exist a single eigenvalue Ei, whose degeneracy, Di, is of 
the order 0{N), whereas all others are of the order 0(1) or 
less. By writing 



1 

N 



Ae-'^'*+ ^ Dn 
one obtains, up to order 0{1/N^) lf74ll : 



-iE„t 



(34) 



iV2 



\di+ J2 D.n2cos[{Ei 



En)t] \. (35) 



The first term on the right-hand side of eq. ( [35T l is of order 
0(1), while the second term is of order 0(1/A^). There- 
fore, for few highly degenerate eigenvalues, the lower bound 
|a(t)|^ will not show a decay to values which fluctuate about 
1/A^ but will rather fluctuate around 1 — 1/A^ at all times. 
Also, 7f(t) will not decay but will fluctuate around the same 
value, since is a lower bound. 

An example of a system with a single highly degenerate 
eigenvalue is the star, having one core node and — 1 nodes 
directly connected to the core but not to each other, see Fig.|2l 
The eigenvalue spectrum has a very simple structure. There 
are 3 distinct eigenvalues, namely Ei = 0, E2 ~ 1, and E3 = 
N, with degeneracies Di ~ 1, D2 = N — 2, and = 1, 
respectively. Therefore, one gets ll74ll 



Pit) - 1 [l + (7V-2)e-* + e- 



-(N-2)t 



1 



1 + (N - 2)e-'^ + e- 



-\(N~2)t 



(36) 



(37) 




FIG. 2: Star with iV - 1 = 12 arms. 
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(b) Star with N=51 



time t 



— 1 1 — I — I I I I I 



T 1 1 — I I I I I I 




— 1 1 1 — I I I I 



1.0 



\ I I mill I I I mill I I I I III! 



J I I I 



_l I I I I I I 



10 



10" 10' 
time t 



FIG. 3: (Colour on-line). Panel (a): p{t) and \a{t)\^ with the 
appropriate scaling t~^^'^ and , respectively, for a ring of size 
iV = 1000. Panel (b): p{t),Tf{t), and \a{t)\^ for a star with iV = 51 
nodes. The inset shows a close-up of 7f(t) and in the same 

time interval. From (74ll . 



Obviously, only the term |(A^ - 2) exp(-it)p/7V2 = (A^ - 
2)'^/N^ in Eq. ([37]) is of order 0(1). AU the other terms are of 
order 0{1/N) or 0(1/A^^) and, therefore, cause only small 
oscillations (fluctuating terms) around or negligible shifts 
(constant terms) from (A^ - 2)^/N^ wl-l/N. 

Figure [3] shows the temporal behavior of 7f(t), |a(t)p, and 
p{t) for a ring (uppper panel) and for a star (lower panel). The 
difference, especially for CTQW, is dramatic. While in both 
cases the CTRW decay to the equipartition value 1 /N, this is 
not so for the CTQW. For the ring \a{t)\'^ (thus also 7f(t)) fol- 
lows the t^^ decay at intermediate times and oscillates around 
the long-time average at long times. For the star, on the other 
hand, there is no decay for both and Tf{t). There are 

only some oscillations around the value 1 — 1 /N. 

In the two examples above the efficiencies of the two pro- 
cesses are vastly different. For the ring there is a considerable 
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FIG. 4: Spidemet graph of generation g = 3 with A*' = 53 = 1 + 
4ELi 3*"' nodes. 



increase in efficiency for CTQW, because for it the probabil- 
ity to return to the origin decays much faster than for CTRW. 
The contrary is true for the star. Here, the CTRW still decay 
while the CTQW remain close to unity. This implies that the 
probability to visit other places than the initial node is - on 
average - very low. Note that 7f(i), |a(i)p, and p{t) are av- 
eraged over all nodes of the network. For the ring all nodes 
are equivalent, but the star has the core node which stands 
out from the rest. Therefore, the temporal evolution is very 
different when the CTQW starts at the core or at one of the 
peripheral nodes. When starting at the core, due to the rota- 
tional symmetry of the network, the walk can be mapped into 
a walk involving only two nodes. When starting at the periph- 
ery, there is no such simple mapping. Since all but one node 
are at the periphery, they lead to the poor performance of the 
CTQW when compared to the CTRW. 

Salimi studied the situation over the semi-regular spider- 
net graph iFttIi . This graph is build in a hierarchical, radially 
symmetric manner, i.e., one starts from a single core node and 
builds up the network generation after generation, as indicated 
in Fig. 13 The procedure used lf73ll74ll77ll also leads to power- 
law behaviors for both CTQW and CTRW fff]. When one 
starts at the core node, one can map the dynamics onto a line, 
where all the states corresponding to the nodes belonging to 
the same generation are summed up to form a new state, rep- 
resenting this generation. Salimi showed that in this case the 
return to the origin (core node) follows a power-law which for 
CTRW goes as t"^/^ while for CTQW it goes as t-^. There- 
fore, a behavior similar to that discussed around Eqs. ( [32] i and 
(l33T l is also found here, since the exponent of the CTQW de- 
cay is twice as large as the exponent of the CTRW decay. 



IV. CTQW ON NETWORKS 
A. Deterministic networks 

1. Two dimensional regular networks 

CTQW on regular 2d structures carry over many of the 
properties of their Id counterparts, since 2d regular networks 
ca be envisaged to be the direct product of two 2d structures, 
vide infra. However, some care is in order, since the symme- 




FIG. 5: Sketch of a square network arranged as a regular lattice with 
the appropriate numbering of the nodes. Note that the actual geomet- 
rical realization can be much more flexible, see text for details. From 

tries observed in Id are not always observed in 2d iItTIi . 

Consider now 2d regular structures of side length N, thus, 
they contain A^-^ nodes giving rise to iV^ basis states IITlll . In 
a pair notation one sets \j) = \jx,jy), where and jy are 
integer labels in the two directions, with jx,jij G [1- see 
FigS This labeling of the states is not to be confused with the 
labeling of the adjacency matrix. Note that capital bold letters 
denote matrices, while small bold letters denote the nodes and 
the states. 

The focus in solid state physics is on systems where Born 
- von Karman periodic boundary conditions (PBC) are as- 
sumed. Now, for an internal site of the network (not on 
an edge or in a corner), the Hamiltonian acting on a state 
b') = \jxJy) reads 

H\ix,iy) = '^\ix,jy) -\ix + l,jy) -\ix-l,iy) 
+ 2|>, jy) - \jx,iy + 1) - \jx,jv - 1>- 

(38) 

PBC extend this equation to all the sites of the network by 
interpreting every integer jx and jy to be taken modulus N . 
With this generalization, the time independent SB 

H\^e)=Eg\^0) (39) 

admits (similar to the Id case) the following Bloch eigenstates 

1*^) = ]^ E (40) 

as solutions, where 9 ■ j stands for the scalar product with 
9 ~ {6x, Oy). The usual Bloch relation can be obtained by 
projecting on the state \ j) such that ^'e^') = O'l^e) = 
e-«(e-j)/iv, thus ^g{jx + l,Jy + 1) = e-*(«-+««)^'e(j,, Jy). 
The PBC restrict the allowed values of 9. In the present case 
ill (side length N), the PBC require that *e(iV + l,jy) = 
^'e(l,jy) and 1'e(j^,7V + 1) = ^-eO':.,!). It follows that 
one must have 6^ — 2mr/N and 9y = 21'k/N, where n 
and I are integers and n^l S [1, -^^]. It is now a simple mat- 
ter to verify that the \^e) also obey {"^eY^e') = 5e e' and 

Furthermore, from Eqs.([39l) and ( |40l i the energy is obtained 

as 

Eg = A-2cosex-2cosey^Ee^^Ee^, (41) 



9 



with Ee^ = 2 — 2 cos 61:1 and Eg^ = 2 — 2 cos 6y. Under 
PBC the two-dimensional eigenvalue problem separates into 
two one-dimensional problems. 

The transition amplitude at time t from state \j) to state \k) 
is now, using Eq.(l40ll twice, [7 1.1 : 



e,e' 

^^e-'^'-'e-^^C^-^') (42) 



In the limit N oo, the sums in Eg. (1421) may be changed to 
integrals; by making use of Eq.dTTli one obtains 



lim akj{t) -— - 



d9x e~*^="('^="^^=")e*^*'^°'^' 



X / d6y Q-i^vi^y-Ov) gi'^i'^osSy ^ 



(43) 



such that 



lim =*'=--^"-z'=='-^«e-'4*J,_,j2i)Jfc„-,„(2t) 

iV — >oo 

(44) 

where J„ (a;) is the Bessel function of the first kind OTSIl . Thus, 
on a network topologically equivalent to a square lattice with 
PBC the transition amplitude between the nodes j and fc is 
given by 



lim nu,j{t) = [Jk.-i.m.h^^,^{2t)f 

iV — >-oo 



(45) 



Now, for finite networks the long-time average Xk.j gives 
more insight into the dynamics for different network sizes N . 
In the 2d case one has 17111 



\k i = lim — dt 



^(fc|e-«*|q„)(gnb-; 



^(fc|q„) (q„ \j) {j\q^m) {qm\k) 



T^oo T J 



(46a) 



^^\^,\^{k\qn){qn\j){j\qni}{qra\k}. (46b) 



One notes that the integral in Eq. (146 al l equals 1 if A„ — Xm 
and otherwise, i.e., it equals (5a„,a,„ ■ Given that some eigen- 
values of H are degenerate, the sum in Eq.( |46b| ) can contain 
terms belonging to different eigenstates |q„) and \qm). Equa- 
tion (I46bl l provides a numerically very efficient way of com- 
puting the Xk,j- Remarkably, one finds that the Xk,j depend 
in an unexpected way on the exact value of the size N of the 
finite network under study. 





FIG. 6: LPs Xk,c to be at node k when starting at the comer node 
c = (1, 1) for networks of sizes (a) iV = 5, (b) iV = 14, (c) iV = 23, 
and (d) N = 47. From HM- 



(a) N=6 




(b) N=15 




FIG. 7: LPs Xk,c to be at node fc when starting at the corner node 
c = (1, 1) for networks of sizes (a) iV = 6, (b) iV = 15, (c) iV = 24, 
and (d) A*' — 48. One may note the asymmetries by comparing to 
Figll From (H. 



When starting at a corner node c = (1,1), one often 
finds that the LPs for the starting node and its "mirror" node 
oc = {N, N) are equal. Figure |6] shows the Xk,c obtained by 
going from the corner node c — (1,1) to the other nodes for 
networks of sizes iV = 5, iV = 14, = 23, and N = 47. 

However, for some particular network sizes the distribu- 
tions of the LPs turn out to be asymmetric. For instance, for a 
network of size N ~ 15 the LP Xoc,c for the CTQW starting 
at node c to be at the opposite corner node oc is less than the 
LP Xc,c to be at the initial node. The same is true for the nodes 
along the edges of the network. Figure Q shows that such 
asymmetries occur for networks of sizes = 6, = 15, 
N = 24, and = 48 (the asymmetries are best seen by look- 
ing at Xc,c and Xoc,c)- The smallest network where asymme- 
tries in the distribution of the LPs are detected has A^ = 6. The 
next ones are found for A^ = 12, 15, 18, 21, 24, 30, 36, • • • . 

An asymmetric LP distribution is particularly evident in the 
difference between Xc,c and Xoc.c- Thus, as an overview 
FiglH] presents as a function of A^ a plot of the (xc,c — 
Xoc,c)N^ values obtained. Note that all N values in Fig|8] 
for which (xc,c — Xoc,c) 7^ are divisible by 3. However, the 
converse is not true, one finds symmetric LP distributions for 
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(a) Rectangle MxN with N=15 
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FIG. 8: Differences between the LPs for CTQW that start at c = 
(1, 1) to be at c, Xc,c, or to be at its "mirror" node oc — {N, N), 
Xoc.c, as a function of the network size A'^, for 1 < < 60. From 



the networks with iV = 3, 9, 27, 33, 39, • • • . 

The above effect is also observed for non-square networks, 
where the number of nodes in the two directions are not equal 
ifrgll . While for periodic boundary conditions in both direc- 
tions there is no asymmetry, there are observable asymmetries 
both for periodic and for open boundary conditions in only 
one direction as well as for open boundary conditions in both 
directions. 

One now fixes iV to be TV = 15 and varies M, taking 
4 < M < 30. Figure |9| a) shows for rectangles the differ- 
ence between the LPs on the initial corner Xc,c and on the 
corner Xoc,c- In the analysed range, 4 < Af < 30, the value 
Xc,c — Xoc,c displays varying patterns. For 4 < M < 13 
one can associate Xc,c ~ Xoc,c = to odd values of M and 
Xc,c — Xoc,c 7^ to even values of M; for larger M the 
situation becomes more complex. Figure|9fb) displays the sit- 
uation for cylinders with = 15 and 4 < Af < 30. Here the 
plot shows Xc,c - Xcy,c, with Cy = {l,N). Figurel^c) shows 
the situation for cylinders with M = 15 and 4 < < 30. 
This last case looks quite regular, with non-vanishing values 
only for A^ = 10, 15, and 30. Hence, Figs. Hb) and^c) 
{N < M) show for cylinders that changes in radius lead to 
more asymmetric situations than changes in length. 

An indication of the origin of the asymmetries is obtained 
by briefly reviewing some of the details of the calculations of 
the LPs. Since the Hamiltonian of the problem separates in the 
two directions, it is easy to show that the transition probabili- 
ties TTk,j{t) can be written as the product of the two separate 
probabilities for each direction 1,79.1 . i.e. 



(b) Cylinder MxN with N = 15 
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(c) Cylinder MxN with M=15 
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FIG. 9: (a) Rectangles MxN with A^ = 15 and varying M: Differ- 
ences between the LPs for CTQW that start at c = (1, 1) to be at c, 
Xc,c, and to be at oc = {M, N), Xoc,c- (b) Cylinders M x A' with 
A' = 15 and varying A/: Differences between the LPs for CTQW 
that start at c = (1, 1) to be at c, and to be at the node Cy — (1, A^). 
(c) Cylinders MxN with M = 15 and varying A^: Differences 
between the LPs for CTQW that start at c = (1, 1) to be at c, and to 



be at the node Cy = (1, A'^). From 1220. 



Now, according to Eq. ( 146 a| ). the LPs are given by I? 



Xk,j = lim - / dtTTk^^jAt) T^ky,jy{t) 

T~>-oo 1 Jq 



1 r 

EFu T lim — / 



dt 



X exp |^-«i(Ae^ - Ae^ + Xq^ - Xe'^] 



(48) 



(47) 



with ak^jAt) = Ee, exp(-iAeJ)(fc^|*0j(\E'eJj:E) , and 
similarly for the y-direction. 



where Fkj is a time independent function, which depends on 



the eigenstates associated with 9^ 



9,,, and 



Because 



of the limit in the time integral in Eq. ( |48T l. there are only 
contributions to Xk,j if a value (Ag^ — Xg' ) for the x-direction 
has a counterpart —{Xg^ — Xg' ) in the y-direction. 

A careful analysis of the differences {Xg^ — Xg' ) indicates 
where the asymmetries stem from. For finite chains one ob- 
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tains (Ae^ — Xq') = 2cos0^ — 2cos6'2.. For simplicity one 
considers now finite N x N networks with OBCs, see iFtTIi . 
because then the eigenvalues are the same in both directions. 
It turns out that for dx ^ d'^ the value (A^^, — Xgi ) appears 
only once or twice for all symmetric cases. However, for the 
asymmetric cases, some of the {\e^ — Xe' ) values (again for 
6x 7^ 0'^) appear more than twice. Therefore, there are more 
contributions to Xkj in the asymmetric cases than in the sym- 
metric cases. 



2. Star-like networks 



(a) 6-fold finite path star graph (b) 4-foid star square graph 




FIG. 10: (a) Star with 6 arms composed of finite segments of length 
4. (b) Star with 4 arms composed of circles of length 4. The dashed 
lines are guides to the eye, connecting all nodes belonging to given 
strata, see text. 



Networks which are not regular but, yet, have symmetry 
properties which can be exploited, are star-like networks. As 
akeady mentioned in Sec. 1111 Al and also discussed in ITSlpTill . 
these networks consist of a central core node (with label 1) 
to which each of the remaining — 1 nodes is attached by 
an individual bond, see Fig.|2] The Hamiltonian has then the 
following structure 



H 



N 

(A^-l)|l)(l|+^(|j)0-|-|l)(j|-b-)(l| 

J=2 



(49) 



Using the Gram-Schmidt orthonormalization procedure, Xu 
obtained the eigenstates of H lls^ 



n+l 



4f j: 



m=2 

for n < A^ 



(50) 



for n ~ N — 1 



1 1 ^ 

\Vn ^ 

V m— 1 



for n — N 



This yields analytic expressions for the transition probabil- 
ities, e.g.. 



7^1,1 W 



N'^-2N + 2 2(A-1) 



iV2 



A2 



cos{Nt). (51) 



One easily verifies that for cos{Nt) = 1, i.e., t = 2TTr/N 
(where r is an integer), there is a perfect revival. Moreover, 
for t — (2r + 1)tt/N all the probabilitiy is distributed over all 
but the core node. This is not true if one starts the CTQW at 
any of the other (A^ — 1) nodes, then one has 



(A"* - AN^ + 5A2 - 2A ^ 

+ (2A^^ - GA^ + 4A) cos(t) 
+ {2N^ - m) cos((Ar - l)t) 

+ (2Ar - 2) cos(At)j 



A2(A^- 1)2 



(52) 



where, without loss of generality, the node 2 was chosen as 
initial node. One sees, especially for large A^, that the prob- 
abilities are mainly localized on the initial node. For all ini- 
tial conditions, the classical transition probabihties pkj (t) ap- 
proach the equipartition value 1 /N. 

Having determined the probabilities to return or to still be 
at the origin allows to calculate the average return probability 
7f(i), see SecHn] One has thus 



N 



N 



fe=l 



1 

N 



ni,,{t) + {N ~ l)TT2.2it) 



which leads to the results given in Eq. dJTl ). 



(53) 



The strong dependence on the initial conditions also carries 
over to the long-time averages Xk.j- One obtains namely 



Xia - (A^ - 2A^ + 2)/A^2^ 
X2a = 

X2,2 = {N'^ - 4:N^ + 5A^^ - 2Af + 2)/ [N'^ (A - 1)^] , 
X3,2 = 2(A2 - a + 1)/ [N^N - If] . 

(54) 



Thus, for large A^ the probability will be concentrated at the 
initial node. 

In a slightly more general setup, Salimi considered CTQW 
on networks (called star graphs in |81]) with are built from 
several sub-networks such that all sub-networks share a single 
node, see Fig.jTO] The adjacency matrix is then a direct prod- 
uct of the separate adjacency matrices with additonal entries 
for the newly created bonds. 

By considering so-called strata, i.e., sets formed by all the 
nodes at the same chemical distance from the central node, 
see the dashed lines in Fig.jTO] Salimi calculated the transition 
probabilities to go from the central node to the different strata. 
For a star with A^ arms each of which having 2 nodes, see 
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Fig. 2 in lISlll . he obtains for the transition probabilities 

2 



^1,1 W = 

and7r3,i(t) = 



1 + iV cos(ViV + It) 



Ar + 1 



iV sin( V^VTTi) 



VN + T 



N{1 - cos{y/N + It) 



{N + 1)VN 



(55) 
(56) 
(57) 



In the case of a star made out of rings of length 4, such that 
all rings share a single node, there are again 2 strata and the 
central node, see Fig. fTOf b) and L81.1 . The transition probabil- 
ities are given by 



and7r3,i(t) = 



l + iVcos(^2(7V + l)t) 



N + 1 



Nsm{y/2{N +\)t) 



^/NTl 



Nil - 2 cos(V^VTTt) 



2{N + 1)^/N 



(58) 
(59) 
(60) 



results very similar to the above. 

For strata more distant from the core, one can still com- 
pute the transition probabilities for different N . Evidently, 
the = 1 case corresponds to a semi-infinite line, while the 
N — 2 case is equivalent to the infinite line, i.e., the transition 
probabilities are given by Bessel functions. 

Letting TV go to infinity, Salimi showed that in all cases 
considered above the transition probabilities reduce to 



7ri,i(t) = cos^t, 
7''2,i(0 = sin^t, 
and7r3,i(t) = 0, 



(61) 
(62) 
(63) 



which is equivalent to the result for a dimer, i.e., a complete 
graph consisting of two nodes. 



probabilities he then obtains 

N^-2N + 2 2{N-1) 



iV2 



7V2 



cos(iVt) 
for k — j 



for k 7^ j. 



(65) 



Therefore, the transition probabilities have exactly the same 
form as the ones for the star graph, when the excitation starts 
at the core. 



4. Dendrimers 

Star graphs of length 1 can also be viewed as being den- 
drimers (Cayley-trees) of first generation. The stucture of 
dendrimers is exemplified in Fig. [TT]for dendrimers of gen- 
erations G = 2 and G — i, with functionality / = 3, see 
also ll82ll . In general, the functionality / gives the number of 
bonds emanating from each node; the generation G refers to 
all the nodes whose shortest distance (in bond units) from the 
central node is not larger than G. Note that the number of 
nodes belonging to the g-th generation (where G > g > 1) 
is 3 • 2^^^ and that it grows exponentially with g. Moreover, 
the total number of nodes in the dendrimer of generation G is 
N ^■i-2° -2. 

The connectivity matrix of these dendrimers has a very sim- 
ple structure. One has An = 3 for all the nodes in the first 
G — 1 generations and An = 1 for the nodes in generation 
G. The bonds are represented by the off-diagonal matrix ele- 
ments Aij . Here, every node in generation g > 1 is connected 
to two consecutively numbered nodes in generation g+1 and 
to one node in generation .g — 1 . 

The eigenmodes of such dendrimers were studied in jssll . 
There, for the dendrimers of generations G = 1 and G = 2, 
the eigenvalues and eigenvectors of A were explicitly calcu- 
lated. The eigenvectors determine the eigenmodes of the den- 
drimer, see e.g. Fig. 3 of jssll . It was further shown that there 
are G + 1 nondegenerate eigenvalues, one of which is always 
Ao = 0. 



3. Complete graph 

The complete graph (for finite N), where all nodes are mu- 
tually connected with each other, shares some properties with 
the star graph. The Hamiltonian now reads 

N 

H={N-l)Y^\^){j\~Y.\^){k\. (64) 

The graph has two distinct eigenvalues Em ~ and En — N 
for n = 1, . . . , N — 1. Xu showed that, similar to the star 
iFtsI Ft?. 80], one can calculate the eigenstates using the Gram- 
Schmidt orthonormalization procedure 1 80]. For the transition 



(a) G=2 (b) G=3 




FIG. 1 1 : Dendrimers of functionality f — 3 and generation G = 2 
(left) and G = 3 (right). From L82J. 
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(a) G=3 
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FIG. 12: Limiting probabilities for dendrimers with (a) G = 3 and 
(b) G — 4. The white lines indicate the limiting distributions for the 
dendrimers of generations smaller than G. From t82i1 . 



When an excitation starts at the central node 1, the dynam- 
ics of this excitation over the dendrimer can be mapped onto 
a line. Remarkably, for the G — 2 dendrimer the transition 
probabilities are fully periodic when the coherent excitation 
starts from the central node 1 (the same holds for the G = 1 
dendrimer, too). Note that due to rotational symmetry, the 
transition probabilities from the central node to nodes belong- 
ing to the same generation are equal. Because of this one 
only has to list three different transition probabilities. It fol- 
lows that there is a perfect revival of the initial state, which 
resembles results obained for continuous ll40ll4Tll and discrete 
quantum carpets |72, ,84,] . 

If the initial excitation is placed at one of the nodes of the 
outermost generation g = G of the dendrimer, the picture 
changes. Classically, the propagation through the dendrimer 
gets to be much slower than in the previous case, see e.g. 
Il58ll85i[86ll . Nevertheless, eventually the excitation will clas- 
sically propagate through the whole graph and in the long time 
limit the probability will be equipartitioned among all nodes. 
Quantum mechanically this effect is even more dramatic. The 
main fraction of TTk,j{t) stays in a small region closely con- 
nected by bonds to the initial node j, and the transfer to other 
sites is highly unlikely. Also at long times the limiting prob- 
ability for the excitation to reach the other branches of the 
dendrimer stays very low. 

Classically, the LP is equipartitioned among all the nodes. 
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FIG. 13: (Color online). Clusters of the same limiting probability 
Xk,j in the branch with the initial excitation for a dendrimer of gen- 
eration G. Nodes connected by thick (red) lines belong to the same 
cluster. From 18211. 



i.e. limt_j.oo Pfc.j (i) = 1/iV for all nodes. Quantum mechan- 
ically this is not the case. Figure [T2l shows the LPs Xk.j as a 
contour plot |82]. Bright shadings correspond to high values 
of the LP, whereas dark shadings correspond to low LPs. The 
diagonal has high values, meaning that an excitation starting 
at node j has a high LP to be found again at node j. The struc- 
tures of the LP distributions of dendrimers are self-similar, 
generation after generation. 

Furthermore, different nodes k and I may have the same LP, 
Xk.j = Xi.j- One hence combines all the nodes having (up to 
our numerical precision, 10~^°) the same LP into a cluster. 
Note, however, that the separation of the nodes into clusters 
depends on the initially excited node, namely on j. For an ex- 
citation starting at the center (node 1), the clusters correspond 
exactly to the different generations of the dendrimer. In the 
general case, when starting from a non-central node, one still 
finds from Fig. [12] that nodes belonging to the same cluster 
also belong to the same generation (the converse is not neces- 
sarily true). 

For larger dendrimers, while the general cluster pattern is 
preserved, some details change. Figure[T3]shows the situation 
for a dendrimer of dimension G = 5, for an excitation starting 
at a peripheral node. Again one indicates clusters by connect- 
ing nodes by thick (red) lines. A change to be noticed is that 
for G > 5 the initially excited node does not form anymore 
a cluster with its next-nearest node of the same generation. It 
appears as if such two nodes only belong to the same cluster 
when the dendrimer has G < 4. Thus the total number of 
clusters is TVc = (G^ + G + 6) /2 for G > 5 and Nc - 1 for 
G < 4. 



J. Husimi-cacti 

As shown in fs?'], it may happen that the excitation occu- 
pies preferentially the bonds between the branching points of 
the dendrimer. Then, the essential underlying structure is dif- 
ferent and is given by sites localized at the mid-points of the 
bonds. The situation is exemplified in Fig.[T4|a), starting from 
a dendrimer of generation 2 (open circles) and indicating the 
mid-points of the bonds by filled circles. Connecting neigh- 
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(a) N=9 



(b) N=21 




FIG. 14: Finite Husimi cacti (filled circles) of size (a) = 9, (b) 
A'' = 21, and (c) TV = 45. (a) also shows with dashed lines and open 
circles the corresponding dendrimer. From Issll . 
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FIG. 15: Limiting probabilites for finite Husimi cacti of sizes (a) 
iV = 21 and (b) N = 45. From LMl- 



boring filled circles by new bonds, one is led to the so-called 
Husimi cactus. Figure [T4l shows three finite Husimi cacti of 
sizes N = 9, 21, and 45. 

In the same way as for the dendrimer, one finds numerically 
that for the finite Husimi cactus consisting of 21 nodes, as 
depicted in Fig. [141 the TPs are nearly periodic when the initial 
excitation is placed on one of the (symmetrically equivalent) 
nodes 1, 2, or 3 of the inner triangle |88]. 

Also in the long time limit, the similarities to the dendrimer 



FIG. 16: (a) Graph consisting of two Cayley trees of generation 
G = 3. (b) horizontal projection of the graph following Ref. 189^], 
(c) vertical projection of the same graph, (d) Vertical projection of a 
similar graph, obtained from two Cayley trees of general generation 
G, indicating the new nodes (clusters) and the dk, see text for details. 
From (23|]. 



are obvious. Figure [15] presents the LPs for two sizes of 
Husimi cacti, = 21 and N = 45. The LP distributions 
are self-similar generation after generation. Furthermore, for 
each size there are LPs having the same value, i.e. Xk.j = Xi,j- 
One collects LPs of the same value into clusters. Depending 
on where the excitation starts, the clustering is different. This 
in analogous to the previous results for dendrimers. Since both 
dendrimers and Husimi cacti lead to similar results, one can 
conclude that here the loops have no significant effect on the 
transition probabilities. 



6. Glued Cayley trees 

One particular example which shares the properties of both 
regular and hyperbranched networks is a network which is 
composed of two Cayley-trees with the same number of gen- 
erations, where the nodes of the last generation are shared by 
both trees, see Fig.[T6l 

Now, depending on the initial condition, the dynamics of 
the CTQW over such a network can change dramatically. 
While the transport from the top to the bottom node is very 
fast 1 89] and comparable to the dynamics on a finite regu- 
lar one-dimensional network, the transport from the left-most 
node to the right-most node is very slow ll23ll compared to the 
transport in dendrimers or Husimi cacti when the excitation 
starts at a peripheral nod e |,82, 88] . 

The authors of Refs. lUland iMl have analyzed CTQW 
over such networks, focussing on walks which start at the 
top node, and looking for the amplitude of being at the bot- 
tom node at time t. The problem can then be simplified by 
considering only states which are totally symmetric superpo- 
sitions of states \k), involving all the nodes k in each row of 
Fig. [161 a), as indicated schematically in Fig. [T6l b). The trans- 
port gets then mapped onto a one-dimensional CTQW l,89il . 
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Other initial conditions for the CTQW are, indeed, possible, 
especially when considering the high symmetry of the under- 
lying graphs. Note that, using for instance the site enumera- 
tion of Fig. [T6l a). a CTQW from node 8 to node 15 is equiv- 
alent to a CTQW from, say, node 10 to node 14. The graph's 
symmetry suggests to collect groups of such nodes into clus- 
ters, while focussing on the transport from left to right. It is 
then natural to view the nodes 8, 9, 10, and 1 1 as belonging to 
the first cluster. The second cluster consists then of the nodes 
4, 5, 16, and 17, all of which are directly connected by one 
bond to the nodes of the first cluster The nodes 2 and 20 of 
the third cluster are all nodes directly connected by one bond 
to the nodes of the second cluster, while at the same time not 
belonging to the first cluster. In general, all the nodes of the 
(fc + l)st cluster are connected by one bond to nodes of the 
fcth cluster and at the same time do not belong to the (fc — l)st 
cluster. 

Let us denote the number of nodes in cluster k by dk- The 
transport occurs now from a cluster to the next, by which the 
original graph gets mapped onto a line in which one new node 
corresponds to a group of original nodes of the graph. For a 
new node at position k G [2, G] one finds that — 2''^'^+^, 
the same being true for the mirror node value, i.e., d^ — 
d2G+2-k- Note that for the end nodes di = d2G+i = 2*^^^, 
the same holds for the nodes next to them. Moreover, for the 
middle node dc+i = 2. 

One now focuses on the transport via the states which are 
totally symmetric, normalized, linear state-combinations for 
all the original nodes in each cluster. Thus, for the fcth cluster, 
whose sites are denoted by n, one has as a new state ll23ll 



G=3 dendrimer 



1 



/dk 



El 



(66) 



The CTQW is now determined by the new Hamiltonian 
H = J A, where the matrix elements of A are obtained from 
the new basis states \ak) and from the matrix A through 



Ajk = {aj\A\ak) 



(67) 



Given the properties of A and the construction of the \ak), 
Eq.(|66l), A is a real and symmetrical tridiagonal matrix, which 
implies a CTQW on a line. The diagonal elements of A are 
given by 



A. 



kk 



ak\A\ak) 

p ^(n'|A|n) = /„ = /fc, (68) 



dk 



nek 



where fk is the functionality of every node in the kth cluster. 
For the sub- and super-diagonal elements of A one finds ll23ll 

Ak,k+i = Ak+i.k = (afe|^|a*:+i) 

= -7rr= E (-'1^1-) 

y'dkdk+l „gfc 

n'ek + 1 

bk 



y/dkd. 



(69) 
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FIG. 17: (a) Transition probability for a CTQW between 

different clusters j of the G = 3 graph. The CTQW starts at the first 
cluster, presented is the situation at times t = 1, 5, 10, 20, 40, 80, 
and 160. (b) Lim iting probability Xji for a CTQW starting at the 
first cluster From 12311 . 



where bk is the number of bonds between the clusters k and 
/c + 1. 

Now, except for the ends and the center of the graph, hk 
equals the maximum of the pair (dfe, dfc+i). Between the cen- 
tral node (dc+i = 2) and its neighbors {do = dQ^2 = 2) the 
numberof bonds is be = 6g+2 = 2. The number of bonds be- 
tween the end node and its neighbor is bi — 69G+1 = 2c?i — 
2^. 

For the graph consisting of 22 original nodes the new ma- 
trix A is a tridiagonal 7x7 matrix, which can be readily di- 
agonalized. The advantage of the procedure is clear: the new 
matrix A depends on the number of clusters and grows with 
(2G+1), whereas the full adjacency matrix. A, grows with the 
total number of nodes in the graph, namely with (3 • 2*^ — 2). 

From Eq.(l66ll the transition amplitude between the state 
jflfc) at time and the state \aj) at time t is given by ll23ll 



ak) 



(70) 



■k+l 



where A is the eigenvalue matrix and Q the matrix con- 
structed from the orthonormalized eigenvectors of the new 
matrix A. 

Now the quantum mechanical transition probabilities are 
given by TTjk{t) — \ajk{t)\'^- Figure [TTta) shows the transi- 
tion probabilities for CTQW over clusters. Remarkably, now 
aheady during short periods of time, such CTQW move from 
one end cluster to the other one. The limiting probability dis- 
tribution, Xjk, which is depicted in Fig. [TTl b), also supports 
this finding. Note that Fig. [TTl b) again reflects the symmetry 
of the original graph. 
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FIG. 18: Dual transformation from Sierpinski Gasket to Dual Sier- 
pinski Gasket of generation g = 3. From 190,1 . 



B. Fractals 

1. Sierpinski Gaskets 

Regular networks have integer dimensions, while some 
tree-like stuctures presented in the previous section, such 
as the dendrimers, can be of "infinite fractal dimension", 
when the (fractal) dimension is taken to be given by ~ 
limi^^oo In N/ In R, N being the number of nodes within 
a sphere of radius R.. In contrast, (deterministic) fractals 
have finite, in general, non-integer dimensions. One partic- 
ular example of a deterministic fractal is the dual Sierpinski 
gasket (DSG) for which the exact spectrum of the eigenval- 
ues of the connectivity matrix is known jg^l . A DSG is an 
exactly-decimable fractal which is directly related, through a 
dual transformation, to the Sierpinski gasket (SG). The DSG 
of generation g can be constructed by replacing each small 
triangle belonging to the SG with a node and by connecting 
such nodes whenever the relevant triangles share a vertex in 
the original gasket (see Fig.fTSl). It is straightforward to verify 
that the number of nodes at any given generation g is = 3^. 

The dual transformation does not conserve the coordina- 
tion number (which decreases from 4 to 3, while the coordi- 
nation number of nodes corresponding to the vertices of the 
gasket remains 2), but it does conserve the fractal dimension 
d f and the spectral dimension d, which are therefore the same 
as for the original Sierpinski gasket, namely df = \n3/ \n2 — 
1.58496... and J= 21n3/ln5 = 1.36521.... The eigenvalue 
spectrum of the DSG connectivity matrix can be determined 
at any generation through the following iterative procedure 
(for more details see |91, 92]): At any given generation g the 
spectrum includes the non-degenerate eigenvalue \n = 0, the 
eigenvalue 3 with degeneracy (3^^^ +3) /2 and the eigenvalue 
5 with degeneracy (3^^^ — l)/2. Moreover, given the eigen- 
value spectrum at generation g—1, then to each non-vanishing 
eigenvalue Xg-i correspond two new eigenvalues accord- 
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FIG. 19: (Color online) Average return probability Ti{t) for the DSG 
of generation g — 5 on a log-log scale. The comparison with the 
classical p{t) evidences that the classical random walk spreads more 
efficiently than its quantum-mechanical counterpart. The dashed line 
represents the envelope of Tt{t). From I.90i1 . 

ing to 
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both A+ and A^ inherit the degeneracy of \g-i- The eigen- 
value spectrum is therefore bounded in [0, 5]. As explained 
in [91], at any generation g, one can calculate the degeneracy 
of each distinct eigenvalue: apart from Xn whose degener- 
acy is 1, there are 2^ distinct eigenvalues, each with degen- 
eracy (3ff-'-i + 3)/2, being r = 0, 1, - 1, and 2'' dis- 
tinct eigenvalues, each with degeneracy (3^^''^^ — 1)/2, being 
r = 0, 1, 5 — 2. As can be easily verified, the degeneracies 
sum up to = 3^. 

For the DSG the CTRW average return probability p{t) is 
readily obtained without numerically diagonalizing the con- 
nectivity matrix, since it only depends on the eigenvalues 
which can be calculated iteratively. Figure [19] displays the 
averaged probabilities p{t), 7f(t) and |a(t)P as a function of 
time, obtained for g — 5. The classical p{t) decays mono- 
tonic ally to the equipartition value 1/A^ , while the quantum- 
mechanical probabilities eventually oscillate around the value 
0.7, which is larger than 3^^. Although the amplitude of fluc- 
tuations exhibited by the lower bound is larger than that of the 
exact value, the agreement between the two quantities is very 
good. In particular, the positions of the extremal points practi- 
cally coincide and the maxima of tt (t) are well reproduced by 
the lower bound. This is analogous to the behavior of walks 
on square lattices, Cayley trees, and stars, as described in the 
previous sections. Notice, however, that for the square lattices 
the lower bound turns out to be exact while for Cayley trees 
and for stars it is only an approximation, which, moreover, 
turns out to be less accurate than what is found for the DSG. 

At short times (t < 57~^) it is possible to construct the 
envelope of 7f(t), which depends algebraically on t 1,90.1 . The 
exponent is sa —0.82, to be possibly compared with d/2 « 
—0.68 which is the exponent expected classically for the in- 
finite DSG. The decay of the average return probability Tr{t) 
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FIG. 20: (Color online) Average return probability vf (t) for the DSG 
of generation g = 4 (bright colour) and g — 5 (dark colour). Its 
lower bound |^ (dashed line) and the classical p{t) (dotted line) 
are also depicted, as shown by the legend. From |9(j. 



for the ST can be estimated as well: its envelope goes like t~'^ 
(classically p{t) ^ t^^, see above), implying a faster delocal- 
ization of the QW over the graph. 

Interestingly, for the DSG, the overall shape of tt (t) does 
not depend significantly on the size of the gasket (see Fig.l20l 
and LSQJ). In fact, the behaviour of Tt{t) is mainly controlled 
by the most highly degenerate eigenvalues. These do not 
change when increasing the fractal size (i.e. its generation). 
These values are: 3 with degeneracy mg(3) — (3^^^ + 3)/2, 
5 with degeneracy TOg (5) = {3^-^-l)/2, (5 ± %/l3)/2 with 
degeneracy 77ig_i(3). 

The inhomogeneity of the pattern of the long-time averages 
Xk.j mirrors the lack of translation invariance of the DSG it- 
self. For instance, being v the label assigned to any vertex 
of the main triangle, Xv.v is a global maximum; off-diagonal 
local maxima coiTespond to couples of connected nodes be- 
longing to different minor triangles of generation g — 1- This 
allows to establish a mapping between the pattern of Xkj and 
the structure of the relevant DSG. 

Since the spectrum of the DSG is known, one can calculate 
the lower bound of the long-time average of the average return 
probability, xib, analytically. At generation g the spectrum of 
the connectivity matrix displays N distinct eigenvalues, where 



N 
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9-2 

J2r 



1 = 3 X 



(72) 



Let us denote the set of distinct eigenvalues by {Xi}^^-^ ^. 
Being m(Ai) the degeneracy of the eigenvalue Aj, one can 
write 



N N 



i=l 



(73) 

Now, one goes over to the space of distinct degeneracies, each 
corresponding to a number p of distinct eigenvalues and gets 
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FIG. 21: Limiting probabilities for the DSG of generation (a) (? = 3 
and (b) g = A, whose volumes are N = 27 and A'^ = 81, respec- 
tively. The white lines enclose the limiting distributions for gaskets 
of smaller generations. Notice that the global maxima lay on the 
main diagonal and correspond to j = 1, 14, 27 and to j = 1, 41, 81 
for g = 3 and for g = 4, respectively. From i90ll . 



the final, explicit formula |19| 
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such that X > 1/3^. Interestingly, in the limit g 
LTA X is finite ll90ll : 
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X > lim Xlb = 77, 



and Xlb reaches this asymptotic value from above. 



(74) 



00, the 



(75) 



18 




<tk|(»)>„ 



FIG. 22: Sketch of a SWN of size AT = 16 containing B = 11 
additional bonds. From [gS]. 



C. Statistical networks 

1. Small-world networks 

The above examples dealt with deterministic networks. 
However, many real systems have stochastic features; thus the 
connectivity can be random. To model this, one can disrupt 
the periodicity of regular patterns by randomly including B 
additional bonds into the network |93]. In such a way one 
creates "shortcuts" and a walker can find shorter paths be- 
tween pairs of sites than on the regular network. So-called 
small-world-networks (SWN) are created by randomly adding 
bonds to a regular one dimensional ring, see Fig.|22] Here we 
do not consider self-connections, i.e., bonds connecting one 
node with itself. 

The general behavior of CTQW on SWN can be analyzed 
by averaging over distinct realizations R 



1 ^ 

-Ti 



(76) 



where the index r specifies the rth realization of the quan- 
tity in question. In so doing one obtains statistical results 
which allow for a comparison with the deterministic situa- 
tion. In particular, we consider here the realization-averaged 
transition probabilities {TTkj{t))R, the averaged probabilities 
(7f(t))fl, their lower bound {a{t))ij, and their classical analog 
{p{t))R. Furthermore, we also discuss the long time average 
(LTA) of each of these quantities: 



dt ■■■ 



(77) 



In the absence of any additional bond, the excitations travel 
along the ring and interfere in a very regular manner, produc- 
ing discrete quantum carpets [72]. Typical for these carpets 
is that they show, depending on TV, full or partial revivals at 
specific times [72]. 

For SWN the situation is quite different. Already a few 
additional bonds obliterate the quantum carpets; the patterns 
fade away [93.1 . By adding more bonds, only the initial node 
retains a significant value for {TTjj{t))ji at all times t. Fur- 
thermore, already for SWN with as little as i3 = 5 the pat- 
tern of {iTjj (t))_R becomes quite regular after a short time, see 
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FIG. 23: Time dependence of the averaged transition probabilities 
(TTfej {t))R for SWN of size iV = 100 with (a) B = 1, (b) B = 2, (c) 
B = 5, and (d) B = 100. The initial node is j — 50 and the number 
of realizations is i? = 500. From I9J1. 



Fig. l23T c). This almost regular shape is reached very quickly 
when B gets to be comparable to N [Fig. l23l d)1. One notes, 
however, that particular realizations may still show (depend- 
ing on their actual additional bonds) strong interference pat- 
terns. These features are washed out by the ensemble average, 
so that only the dependence on the initial node stands out. 

Since CTQW on SWN always carry the information of their 
initial node j, the averaged probabilities to return to j are a 
good measure to quantify the efficiency of the transport on 
such networks, see iItsIi . 

Figure |24] shows in double-logarithmic scales the ensem- 
ble averages {p{t))R, (tT (t))^,, and {a{t))ii for SWN with 
N ~ 100 nodes and 5 = 1,2, 5, and 100. For classical trans- 
port [Fig.|24ta)] the initial decay of {p{t))]i occurs faster for 
larger B. The decay at intermediate times follows a power-law 
(i~^/^) for the ring (as is clear from the linear behavior in the 
scales of the figure) and changes to a stretched exponential- 
type when B is large [94,1 . Thus, a classical excitation will 
quickly explore the whole SWN, so that it will occupy each 
site with equal probability of 1/N already after a relatively 
short time, see the final plateau in Fig.|24ta). 

Quantum mechanically, however, the situation is more 
complex. Fig. |24jb) shows the ensemble average (tt (t)) r. 
For a ring of N nodes and for times shorter than roughly 
N/2 {TT{t))]i displays a quasiperiodic pattern (black curve), 
the maxima of which decay as t^^. At longer times interfer- 
ence sets in and leads to an irregular behavior at times longer 
than N/2 [73]. Now, for SWN, as long as B is consider- 
ably less than N, the periodic pattern still remains visible; in 
Fig. \24\ h) one can follow how an increase in B (red, green, 
and blue curves) is smoothing out the curves, so that both the 
heights of the first maxima and the depths of the minima de- 
crease. At longer times the SWN patterns are flattened out 
and {Tf{t))R tends towards a limiting value. With increasing 
B this asymptotic domain is reached more quickly, such that 
for larger B the crossover from the quasiperiodic behavior at 
short times to a smoothed out pattern at longer times is shifted 
to smaller t. 

Figure \24\c) shows the lower bound of 7f{t), namely 
{\a{t)\'^)ii averaged over the realizations. One notices that 
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FIG. 24: (Color online) Time dependence of the averaged probabil- 
ities (a) {p{t))R, (b) (7f(t))_R, and (c) {\a{t)f)R for SWN of size 
N = 100 with B = 1,2, 5, and 100. The number of realizations is 
R = 500. From 191 



the overall behavior of Figs.l24Tb) and l24Tc) is quite similar. 
However, the limiting values at long times differ For the LTA 
of {W{t))ii one has (see also Eq. (17) of Ref. [71,1 ) 
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of course, not valid for the ring, see below). In Eq. (|79] l the 
triple sum adds then to RN, so that the rhs equals 1/A^. On 
the other hand, Eq. ^ leads to Q 



{X)i 



1 

RN 



(80) 



This expression depends on the eigenstates; in fact the rhs of 
Eq. dSOl ) is the ensemble average of the average participation 
ratio of the eigenstates |<i>„.r)- Equation dSOl l is well known 
in the theory of quantum localization, see, e.g.. Sec. V. A. in 
ll95ll . For the ring the eigenstates are Bloch states. 
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from which |(fc|<i>„)| = follows for all |(f>„). 

Now, increasing 13 results in an increase of {x)r, starting 
from the corresponding value for the ring {B = 0, only one 
realization, and N even) 
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where Xn = (2^^ - 2)/N'^. Equation dZDl yields a 1/iV de- 
pendence for the LTA of (|a(t)p)fl, which by rescaling with 
(Xri"<r)fl: ^ 1/^ would result in a constant value for large 
N 1 93]. However, rescaling {'x)r with (Xring)fl shows an 
increase with N of (x) R / iXiin'y) R which is less than linear, 
thus, (x)ij depends on N as l/N", with i/ € [1,2]. 

The fact that (%) r for SWN increases with increasing B 

points towards a change of |(fc|<I>„)| from the value l/iV^. 
The situation may be visualized as follows: For the ring all 
eigenstates are Bloch states and hence are completely delo- 
calized. Going over to SWN and increasing the number of 
additional bonds B leads to localized states at the band edges 
and to fairly delocaUzed states well inside the band. The in- 
crease of (x) R. is thus mainly due to the localized band edge 
states. 
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where 5{En.r — En>,r) — 1 for En,r — En',r and 6{En,r — 

En' .r) ~ otherwise. For the long-time values 

for different B collapse to one value. In fact, the LTA of 

(|a(t)|2)^ obeys 
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(79) 



as can be immediately inferred from Eq. ( l24l l. Thus this quan- 
tity is only a function of the eigenvalues En.r and does not 
depend on the eigenstates |$n.r)- In order to quantify the dif- 
ferences between Eqs. ( |78] l and ( |79] l for SWN, one assumes 
that all the eigenvalues are nondegenerate (this assumption is. 



2. Erdds-Renyi networks 

Somewhat similar to the SWN is the Erdos-Renyi network 
(ERN). One starts with N disconnected nodes, every pair of 
nodes is then connected with the probability p, where only 
single connections between two nodes are allowed. One can 
in turn associate to p an average degree k of the nodes, which 
is related to p by fc = p{N — 1). For large network sizes N, 
the degree distribution P{k) of the ERN is Poissonian peaked 
at 1. 

Xu and Liu showed in Ref. f9^ that in the ensemble aver- 
age, the average return probabilities displtw a behavior very 
similar to the SWN, see Sec. IIV C l1 and ll93ll . Although the 
probabilities decay, they do so only until they reach a plateau 
(for N = 100), at a level considerably higher than the equipar- 
tition value of 1 /N. The behavior is only weakly affected by 
the value of the average degree k. 
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FIG. 26: Example of an apollonian network of generation G = 2. 



FIG. 25: Example of a scale-free network with = 36 nodes and 

(t> = 5/2. 



The height of the plateau is determined by the long-time 
behavior of the transition probabilities. Due to the ensemble 
average, all structure of the non-diagonal elements of (^kj) 
disappears and only the main diagonal (S^ j) remains. 



3. Scale-free networks 

The distribution P{k) for the number k of bonds emanating 
from a node does not need to be Poissonian in general. Net- 
works for which the distribution P{k) follows a power-law, 
i.e., P{k) ^ k^'^, are called scale-free networks (SFN), see 
Fig.|25] These have been proven useful in various fields of re- 
search from biology to social sciences. Xu and Liu have con- 
sidered CTQW over such structures |97]. They distinguish 
between deterministic scale-free networks (DSFN) and ran- 
dom scale-free networks (RSFN): the latter may obey distinct 
building procedures which can lead to either tree-like struc- 
tures or structures containing loops. 

For the DSFN, Xu and Liu determined the return probabli- 
ties Trj.j{t) and long-time averages Xkj- They find that there 
is a strong dependence of Trj.j{t) on the initial node j, which 
can even result in (almost) complete revivals of the initial con- 
dition. The transition probabilities translate directly into the 
long-time averages Xk.j, where one can identify clusters of 
nodes having the same Xkj- However, the patterns obtained 
are quite distinct from the previously found patterns for the 
dendrimers ll82ll or for the Husimi cacti ifssll . As mentioned 
earlier, this is a direct consequence of the fact that the Xkj 
mirror the topology of the network. 

For RSFN, Xu and Liu calculated the ensemble average 
over many realizations of individual RSFN (9^ . In the ensem- 
ble average the return probabilities do not oscillate but rather 
reach a stationary value, which differs for different initial 
nodes. In all cases, however, this stationary value is roughly 
one order of magnitude larger than the classical equipartition 
value 1/iV obtained for CTRW. Moreover, the higher the sym- 
metry of the initial node (where on average the central nodes 



have the largest symmetry), the larger is the stationary value. 
Large average return probablities are also found in the long- 
time averages Xk.j, where in the ensemble average the values 
on the diagonal (for k = j) are much larger than the values 
for k^j. 



4. Apollonian networks 

So-called Apollonian networks are models for networks 
which have small-world as well as scale-free properties. 
CTQW on such structures have been investigated by Xu et 
al. [98]. The network is generated from a triangle (N = 3 at 
generation g ~ 0). In generation g — 1 a single node is added 
which is connected with all three nodes of g — 0, this divides 
the triangle of g — into three distinct triangles. In g = 2 
three new noded are added which are placed at the center of 
each of the three triangles of g = 1, see Fig.|26l These three 
nodes divide each triangle of g — 1 into three new triangles, 
such that there are now nine new triangles in total. The iter- 
ation proceeds in the same manner, such that in generation G 
the total number of nodes is iV = 3 + (3*^ - l)/2. 

For G = 1 the Apollonian network is identical to the com- 
plete graph of size = 4, therefore the transition probabil- 
ities for CTQW are fully periodic. It also turns out that for 
G = 2 (N = 7), see Fig. |26] and when choosing as the ini- 
tial node the central, most symmetric node 4, the transition 
probabilities are also fully periodic, being namely ll97ll 



7I'/c,4(i) 



_ f [37+ 12cos(7t)]/49 forfc = 4 
~ 1 [2 - 2 cos(7i)] /49 for k ^ 4. 



(83) 



However, also here Xu et al. observe a strong dependence 
on the initial condition. If the initial node is not a central node 
being, say, node 1, the transition probabilities tti are still 
periodic, although there is no perfect revival. Enlarging the 
networks lets the strong dependence on the initial condition 
be more pronounced. But, similar to the dendrimers ll82ll . Xu 
et al. also identify clusters of nodes which, in the long-time 
average, have the same value for their probabilities Xk,j- 
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V. EXTENSIONS 
A. Systems with long-range interactions 

So far the Hamiltonian for CTQW and the transfer ma- 
trix for CTRW have been directly related to the connectiv- 
ity matrix. However, this matrix is purely topological (in- 
dicating whether or not two nodes are connected) and does 
not take metric (i.e., distance dependent) aspects into account. 
Some adjustment is thus necessary if one considers distance- 
dependent interactions. As an example, consider the net- 
work being embedded in the d-dimensional space; then to ev- 
ery node j one associates a vector j whose coordinates are 
{xji, . . . , Xjd}- The distance between two nodes j and k is 
then given by the euclidean norm R = \k — j\. Consider 
now systems in which the interaction between the two nodes 
j and k decreases with increasing distance. An example for 
such a system could be a collection of dipoles, interacting via 
dipole-dipole forces, whose potential decreases, to a good ap- 
proximation, as R~^. Now, in terms of the connectivity ma- 
trix, a node j is not only connected to its nearest neighbors 
but also to other nodes. Thus, due to the decaying interaction 
potential the transition rates are not the same for all bonds but 
they depend on the distance. 

Take as a first example a one-dimensional network with pe- 
riodic boundary conditions (i.e., a discrete ring) ll99ll . Here, 
when the interactions go as R^'^, the Hamiltonian has the fol- 
lowing structure: 



R- 



(2\n){n\-\n-R){n\-\n + R){n\y 

n=l R=l 

(84) 

where i?inax is a cut-off for finite systems. Note, that in the 
infinite system limit one first takes — > oo before taking also 
^max — >■ oo. For the cases considered here, namely 7 > 2 and 
A^ of the order of a few hundred nodes, a resonable cut-off is 
^^max — N/2, which is also the largest distance between two 
nodes on the discrete ring. In this way, to each pair of sites a 
single (minimal) distance and a unique interaction is assigned. 

For all 7, the eigenstates are again the Bloch states j^'g) 
given above. The fact that the eigenstates are not affected by 
the long-range interactions is due to the translational invari- 
ance along the ring. For other systems without such an invari- 
ance the eigenstates will also change depending on the type of 
the interaction. From the eigenstates one obtains th e eig enval- 
ues which now do depend on the interaction range ||99I1 : 



E-riO) = J2 fi~^[2-2cos(6li?)]. 



(85) 



R=l 



The DOS (E) is obtained by inverting Eq. dSST l and tak- 
ing the derivative with respect to E^. In the NN-case (7 — 00) 
one gets the known DOS 



Poo{E) = (tt^/AE-E^) \ (86) 
For 7 = 2 one can approximate the sum by letting i?,nax 



00, which yields £'2(6') = tt6 ~ ^^/2, and one obtains 

, -1 



P2{E) = {nV2y^Tr^/2- E 



(87) 



In the intermediate range there is an analytic solution for 7 — 

4,namely£:4(6') ^ 9'^ /24:- 7:9^ /6 + T:^e^ /6 (seeEq. 1.443.6 
of |75]), which yields |1000 



P4{E) = 2TT{2/3Y/^JE{nyVM) - E^ 



(88) 



One assumes the following general form for the DOS in order 
to interpolate between P2{E) and p^{E) ll99tl : 



C-yE"' - E^ 



(89) 



with a E [0, 1] and (3 E [1, 2]; is a constant related to the 
maximal energy, = (-E'y,max)'^ 

Having the DOS at hand, the integrals in Eqs. dSOl l and (ISTT l 
can be calculated - at least asymptotically - for large t. In the 
classical case Eq. ( l30b will be dominated by small values of 
E when t becomes large. The DOS yields 



fora = l 



(90) 



Quantum mechanically one knows that for the NN-case 
T^ca{t) ^ see for instance ||73[]. Considering now the 
other limiting case, 7 = 2, one has 19911 



7f2(t) 



dE 



exp{~iEt) 



7rV2v/7r2/2--£; 



(91) 



Thus, the behavior for long times is the same for Tf2{t) and 
7Foo(i)' which suggests that for all one-dimensional lattices 
with extensive (7 > 2) interactions the long time dynamics of 
the excitations is similar, no matter how long- or short-range 
the step lengths are. This is in contrast to the classical case, 
where only CTRW with 7 > 3 belong to the same universality 
class. 

These results are corroborated by analytically evaluating 
TT-y{t) using the stationary phase approximation (SPA) lIlOlll . 
For large A^, oi-y {t) can be written in integral form 



1 



dd eyiY>{iEy{e)t). 



(92) 



The SPA asserts now that the main contribution to this in- 
tegral comes from those points where Ej{9) is stationary 
[dE^{e)/de = E'^(Q) =0]. For 7 = 2, E2(B) has only 
one stationary point in e [0, 27r[, namely ~ tt, leading to 



7f2(t) = |a2(i)|^ 



1 



2^t\E'i(j:)\ 



(93) 



which does not show any oscillations and coincides with the 
long time limit of Eq. ( |9T| i. For 7 > 2, E^{6) has two sta- 
tionary points in the interval 9 E [0,2tt[, namely 9q — and 
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FIG. 27: (Color online) (a) Classical p^{t) and (b) quantum mechan- 
ical TT^ (t) for a discrete ring with TV = 10000 nodes with 7 = 2,3, 
4, and oo. From | 99]. 



It turns out that the Bloch states are a lso th e eigenstates of this 
Hamiltonian: their eigenvalues read Il02ll 



En = 2m — 2 cos(j0„). 



(97) 



Inserting these values into Eqs. JTOt and (|7]) allows to study 
the dependence of CTQW and of CTRW on to. As has been 
shown by Xu ll02ll . increasing m results (as intuitively ex- 
pected) in a faster transport, both for CTQW and CTRW. 
Morevoer, when considering the long-time average Xkj, Xu 
finds characteristic peaks which depend on m. Especially for 
even N, the probability to be at the initial node and the proba- 
bility to be at the exactly opposite node, i.e., the node j±N/2, 
are not necessary equal 1 102], as it is the case for only nearest- 
neigh bor co uplings. Depending on to, these two values may 
differ il02ll : an explanation of this effect is still lacking. 



6*0 = TT. Then a-^ (i) is approximately given by the sum of the 
contributions of the two stationary points. Consequently L99tl . 



1 



1 



1 



2cos{t[E^{0) - E^{tt)] +tt/2} 

y^|i?;'(o)£;;'WI 



(94) 



The classical and quantum mean square displacements 
(MSD) are in line with these findings. Now, the MSD for 
CTRW/CTQW on the discrete ring with initial site j are given 
by||99[ 



(Kit)) 



cl; qm 



1 ^ 



(95) 



fc=i 



where T'gCi) = for CTRW and T'glt) - 7rg(t) 

for CTQW. Figure [27] shows numerical calculations of p^{t) 
and Tf-f{t) for different 7 and a discrete ring of = 10000 
nodes. Clearly, p^^{t) changes when increasing the step width 
fromNN steps to long-range steps, see Fig.lZTl'a). While 
for 7 > 3 decays as t~^^^, the power law changes to for 
7 = 2. In contrast, the decay of the maxima of the quan- 
tum return probability Tf^{t) follows t^^ for all 7, Fig. lZTl' c). 
Long-range steps lead only to a damping of the oscillations 
and to an earlier interference once the excitation has propa- 
gated around half of the ring. 

As an extension of the ring topology just discussed fg^, 
Xu considered networks where not only the nearest neighbors 
are connected but in which each node is connected to its 2to 
nearest neighbors 1 102]. This model differs from the one with 
long-range interactions, where the interaction decreases with 
the distance |A: — j | between the nodes k and j. 

Now the action of the Hamiltonian on state \ j) reads 



H\j) = i2m + l)\j)- 



[1 + z) 



(96) 



B. Systems with disorder and localization 

In real physical systems, under the influence of the sur- 
roundings, the couplings between the nodes may differ In 
a static picture, one can introduce disorder by adding to the 
unperturbed Hamiltonian i?o a disorder operator A, i.e., by 
setting = i?o + A 1110311 . The disorder matrix A = (A/j ) 
is taken to have non-zero entries only at the positions for 
which Hi j 7^ 0. For different strengths of disorder, the el- 
ements A; J ~ Aj^i are chosen randomly (drawn from a nor- 
mal distribution with the zero mean and unit variance, and 
then multiplied by a factor of A which takes values from the 
interval [0, 1/2]). Note that under these assumptions for the 
(static) disorder the connectivity of the graph is essentially 
unchanged, i.e., there are no new connections created nor are 
existing connections destroyed. Therefore, the only non-zero 
matrix elements of H are those of the initial A. The action of 
the new Hamiltonian on a state |j) reads then 



H\3) = 
= 2|j) - 



\3 - 1) - \J + 1) 



A,j-ib-l) 



(98) 



In the following two cases of disorder are considered: 

(A) Diagonal disorder (DD), where A^ ^ 7^ and A/ j = 
for I ^ j- Here, a random number is assigned to each Ajj, a 
procedure which leads to N random numbers. 

(B) Diagonal and off-diagonal disorder (DOD), where a 
random number is chosen for each Ajj and for each Aj 
For this, 2N random numbers are needed. 

Introducing disorder into the system in this way has conse- 
quences for the relation between the CTQW and the CTRW. 
In CTRW the transition rates, given by the entries of the the 
transfer matrix T, are correlated, i.e., for each site the sum of 
the non-diagonal rates for transmission from it and the diag- 
onal rate of leaving it are the same. In the cases considered 
here, a direct identification of the Hamiltonian H with a clas- 
sical transfer matrix is not possible anymore. However, the 
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FIG. 28: (Color online) Long-time ave rage Xkj for different A and 
(a) TV = 100 and (b) TV = 101. From liol . 



DOD and DD Hamiltonians are widely used in quantum me- 
chanical nearest-neighbor hopping models, to which also the 
CTQW belong. Furthermore, we still consider transport pro- 
cesses on graphs which have the connectivity matrix A, but 
the direct connection between H and T is lost. 

Now, consider again rings of TV nodes, where at time t = 
the excitation is assumed to be localized at node j. The 
above system is similar to the Anderson model lITlll . which 
has been found to show (strong) localization around the initial 
condition. The same happens here. In the ensemble average 
as well as for single realizations, an exc itation starting at j 
remains localized in the vicinity of j ll03ll . This effect is best 
seen in the long-time average. Figure l28l shows Xk.j for two 
different ring sizes, N = 100 and TV = 101, and varying 
A. Clearly, the larger is A the more is the long time average 
localized around the initial node j — 50. 

Quantum walks in random linear environments were also 
studied by Yin et al. fl04ll . They calculated numerically 
the quantum carpet structures for DD. It turns out that in the 
course of time the excitation stays localized around its initial 
node; also depending on the strength of the disorder, some 
structure of the original quantum carpets may still remain vis- 
ible. 

Yin et al. also studied dynamic diagonal disorder, where the 
diagonal elements of H are rapidly varying with time and the 
update of {H)jj is done at times comparable or much smaller 
than the time step of the numerically determined dynamical 
changes il04ll . For dynamic disorder, the interference patterns 
making up the quantum carpets are washed out. However, dis- 
tinct from the case of static disorder, no localization can be 
seen. Moreover, after the temporal range in which interfer- 
ence is lost, the dynamics becomes classical. This crossover 
behavior is also observed in the mean square displacement, 
which changes in a certain time range - depending on the 
strength of disorder - from the ballistic behavior ^ P to the 



diffusive behavior ~ 1 1110411 . 



VI. SYSTEMS WITH ABSORPTION 

In general, an excitation does not stay forever in the sys- 
tem in which it was created; the excitation either decays (ra- 
diatively or by exciton recombination) or, e.g., in the case of 
biological light-harvesting systems, it gets absorbed at the re- 
action center, where it is transformed into chemical energy. In 
such cases, the total probability to find the excitation within 
the network is not conserved. Such loss processes can be mod- 
elled phenomenologically by changing the transfer matrix or 
the Hamiltonian |50, 105-107]. To fix the ideas, we consider 
networks in which the excitation can only vanish at certain 
nodes. These nodes will be called trap-nodes or traps. In the 
absence of traps, let the transfer matrix and the Hamiltonian of 
the corresponding network be To and Hq, respectively. Take 
now M out to the TV total nodes to be traps and denote them 
by m, so that m e Ai, with Ai C {1, . . . , TV}. The trapping 
process is now modelled by introducing a trapping matrix F 
which is given by a sum over all trap nodes; T has only diag- 
onal elements, i.e.. 



(99) 



(in the following one assumes that r„i = F > for all m). 

CTRW with decreasing exciton probabilities due to trap- 
ping are well described through the following transfer matrix: 

T EE To - r. (100) 

The total Hamiltonian H corresponding to trapping is then: 

H=Ho iT. (101) 

Note that the connection between CTRW and CTQW is now 
less direct than before. For CTRW the term corresponding to 
trapping has only real elements and the total transfer matrix 
stays real. For CTQW, however, the trapping term has purely 
imaginary elements. As a result, H is non-hermitian and has 
TV complex eigenvalues, Ei ~ ei — i'^i {I = 1, . . . , TV). In gen- 
eral, H has TV left and TV right eigenstates and ($/|, re- 
spectively. It turns out that in most cases the eigenstates of H 
form a complete and biorthonormal set, see, e.g., Ref. [108.1 . 



N 



1=1 



and 



(102) 



Both, eigenvalues and eigenstates, will be different for CTRW 
and CTQW, because the incorporation of the trapping process 
is different. 

If the trapping strength F is small compared to the cou- 
plings between neighboring nodes, perturbation theory allows 
to relate the real part of the eigenvalues to the eigenvalues of 
the unperturbed Hamiltonian Hq. Let be the ^th eigen- 

state and e K. be the kh eigenvalue of the unperturbed 
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system with Hamiltonian Hq. Up to first -orde r, the eigenval- 
ues of the perturbed system are given by lIlOSll 



Ei=E, 



(0) 



(103) 



Therefore, the correction term determines the imaginary parts 
7;, while the unperturbed eigenvalues are the real parts e; — 
Eji^\ Moreover, the imaginary parts are solely determined 
by the contribution of the eigenstates of the network without 
traps at the ti-ap nodes m. 



A. Average survival probability 

In an ideal experiment one would excite exactly one node, 
say j ^ A4, and read out the outcome TTkj{t), i.e., the proba- 
bility to be at node fc ^ at time t. However, it is easier to 
keep track of the total outcome at all nodes k ^ Ai, namely of 
Sfc^A^ '^kj{t)- Since the states |fc) form a complete, orthonor- 

mal basis set one has J^k^M I^X^I = 1 - J2meM \m){m\, 
which leads to [50]: 

N N 



1=1 



LV = 1 

E mv){^v\m){m\^i){^i\j). 



(104) 



By averaging over all j ^ A^, the mean survival probability 
is given by [l50ll 



E E ^^■^■(*) 




(105) 



1-2 E ($/k>(m|$/) 



$r|m)(m|$/) 



(106) 



For long t and small M/N, Eq. ( 11061 ) simplifies con- 
siderably: At long times the oscillating term on the 
right hand side drops out and for small M/N one has 

2Z]meAi(^'l™)("^l*i) < 1- Thus, UM(t) is mainly a sum 
of exponentially decaying terms l.50il : 



1 



AT 



N - M 



1=1 



Asymptotically, Eq. ( 1107b is dominated by the 7/ values 
closest to zero. If the smallest one, 7min, is well separated 
from the other values, one is led for t ^ l/7niin to the expo- 
nential decay found in earlier works, Il]\,j{t) = exp( — 27niint) 

Coi. 



Such long times are not of much experimental relevance 
(see also below), since most measurements highlight shorter 
times, at which many 7/ contribute. In the corresponding en- 
ergy range the 7/ often scale, so that in a large I range one finds 
7/ ^ al^. The prefactor a depends only on F and N |109|]. 
For densely distributed 7; and at intermediate times one has, 
from Eq. ( fTOTt , ||50I1 



nA/(t) 



dx • 



dy 



(2ai)-i//' 



(108) 



Analogously, the mean survival probability for CTRW is 
given by 



PM{t) 



1 



N - M 



E E p^i'^^)- 



(109) 



If the smallest eigenvalue, Ai, is well separated from the rest, 
PM{t) turns very quickly into a simple exponential decay. 
Then, for not too small times, it can be shown that l,105i] 



PM{t) 



1 



N -M 



k^M 



(110) 



B. Regular networks 

1. Ring with traps 

Some of the nodes of a ring of N nodes are now taken to 
be traps. Depending on the particular choice of trap arrange- 
ments, t he ave rage survival probability 11 ^ (t) shows different 
features iflOell . 

As already discussed, the eigenstates of the ring without 
traps are Bloch states. Furthermore, if the trapping strength is 
small compared to the interaction strength between the nodes, 
the imaginary parts of the Hamiltonian of a network with traps 
are given by the contributions of the eigenstates at the trap po- 
sitions, see Eq. ( 11031 ). However, since all eigenvalues except 
El = (and for even TV also except Ej^f2 — 4) are two-fold 
degenerate, some care is in order when applying perturbation 
theory. 

Take now N to be even. Then, for I = 1 and for / = N/2 
one has 



E^ 



(1) 



(m|$ 



E] cxp[— 27ji]. (107) Furthermore, 



El 



iT- 



N 



and E 



N/2 



-iT- 



N' 



For / different from 1 and from N/2 one sets 



(0)\ 



(111) 



(112) 



(113) 
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and applies t he ex pression valid for two-fold degenerate solu- 
tions of ffo iflOell : 



4 



{Vi,i-Vn-ln-i) +m.N-i\' 



1/2 



(114) 



where one takes the positive sign for / G [1, N/2 — 1] and the 
negative sign for I e [N/2 + 1,N — 1]. Now one has 



Vii = Vn~i.n-i 



N' 



(115) 



independently of the trap arrangement and 



M 



Vi^N-i = exp{2^7rm, [I - {N - l)]/N} 



M 



N 



e:cp{4iTrlmj / N) . 



(116) 



Inserting the last results into Eq. ( II 141 ) yields (TOj 



(1) 



M± 



A/ 

^2iTT2lmj/N 



(117) 



Notice that for special trap arrangements the £"1^' can be cal- 
culated exactly: The most striking results are obtained when 
the exponential in the sum in Eq. ( 1117b equals one of the val- 
ues from the set {1, i, —1, —i}. Then the absolute value of the 
sum reduces to | J2j=i ex-p{iATTlmj /N)\ — M. 

Now, for a single trap, which without loss of generality is 
placed at position mj = j = 1, one has e[^^ = —iV/N[l ± 

exp(i47rl/iV)]. Therefore, one has e[^^ = for ? = TV and 
for I = N/2. As a consequence, IlAf (t) will not decay to zero 
but to a constant value given by 1/ [N — 1). 

Another example is a periodic distribution of traps with 
■nij = jN/M, while N/M G N, see Fig. I^a). Then it is 
straightforward to show that if 21 /M E N the sum in Eq. dl 17t 
equals one of the values from the s et |1 , i, —1, —i}. The total 
number of such values is given by [il06ll 



L(iV-2)/AfJ 
[(A^- 2)/2AfJ 



for even M, 
for odd M, 



(118) 



where [a;J denotes the largest integer less than or equal to x. 
In particular, for both M = 1 and Af = 2, |T| = A^/2 - 1. 
Hence, for large structures with M <C N, IlM{t) decays 
asymptotically to 1 /M (even case) and to 1/(2M) (odd case). 
Figure [30] shows results obtained for a ring of size A^ = 300 
with a periodic arrangement of Af — 10 (|T| = 29) and 
M = 75 (|T| = 1) traps. Consequently, the survival prob- 
ability nj\/(t) decays to the constant values 1/10 and 1/225, 
respectively. From a physical point of view, the finite limit for 
the survival probability stems from the existence of station- 
ary states to which the nodes in A4 do not contribute, so that 



(a) periodic arrangement of traps 




(b) sequential arrangement of traps 




FIG. 29: Periodic (a) and sequential (b) arrangements of traps on a 
ring. 
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FIG. 30: Survival probabilities TlM(t) (continuous line) and Pi\/(t) 
(dotted lines) on a ring of size A'' — 300 and F — 0.01 in the pres- 
ence of M = 10 and of M = 75 traps arranged per iodic ally, i.e. 
rrij — jN/M. Note the semilogarithmic scales. From Il06ll . 



they never "see" the traps. This genuine quantum-mechanical 
effect has no counterpart in the classical world where, for fi- 
nite structures, the survival probability always decays to zero 
in the presence of traps. In particular, as shown in Fig. [SO] 
Pm(^) decays exponentially, as expected. 



For a sequential arrangement of traps, such that rrij 
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FIG. 31: Survival probability Da/ (f) on rings of size TV = 32, 48, 64 
and 96 with a sequential arrangement of M = N/2 traps for T = 
0.04, 0.01. 0.004, as indicated.The straight lines represent Eq. l ll21t . 
From ilOai . 



and i = 1, M, see Fig.|29lb). Eq. ll 16l can be written as 



M 



VlN-l 



^ -ir/A^^exp(4i7r/j7iV) 
-iT exp{Ai7TlM/N) - 1 



(119) 



TV exp(47ri//A^) - 1 
-iT suv{2itMI/N) 



cxp{4nil/N) 



N sin(27r//A^) 



exp[2i7r;(M + 1)/A^], 



which yields 



^ — [^^^ sM2nl/N) )■ ^^^^^ 

Notice that since I ^ N/2 and I =^ N then 21 /N ^ N, while 
for 21M/N e N then E^^^ = E^^li = -iTM/N. In particu- 
lar, when M = N/2, ji = M /N for each value of I € [1 , N] . 
As a result, and by neglecting oscillations, one has lll06ll 



M 



^-2TtM/N 



-Tt 



(121) 



N - M 

which is independent of N . As shown in Fig. [31] the expo- 
nential behaviour predicted by Eq. (1121b holds also for inter- 
mediate times. 



2. Line with traps 

An example of a network with traps which allows to study 
different time scales is a finite line of nodes with traps at each 
end. The Hamiltonian is thus li5Qil 

N-l 

H= ^ (^2|n)(n| - |n- l)(n| - |n+ 

11=2 

+ |l)(l|-|2)(l| + |TV)(TV|-|iV-l)(TV| 
+zr(|l)(l| + |TV)(TV|). (122) 
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FIG. 32: (Color online) Imaginary parts 7; (dots) in ascending order 
for = 100 and F = 1. Note the shortened y axis. The inset shows 
7z in log-log scale for Z = 10, . . . , 90. From [50.1 . 



Without loss of generality, an eigenstate of t he fin ite chain 
without traps can be written as (/ = 1, . . . , N) lIlOSll 



I* 



(0)\ 



N 

2 ^ 

77 cos 

N ^ 



for « = iV 



[(2j - \)ei/2\ \j) else, 



(123) 

where for convenience one takes 61 = Tr{N — l)/N e [0, 7r[; 
the corresponding eigenvalues are Ef''' = 2 — 2 cos 9i (note 

that the smallest eigenvalue is = 0). Thus, first order 
perturbation theory yields from Eqs. (I103l l and (11231) as imag- 
inary parts = 2T/N and 7/ = (4r/A^) cos^ [61/2) for 
I = I, . . . ,N - I. Indeed, for / < this means that 7; P, 
so that the average survival probability scales in the corre- 
sponding time interval as TlM{t) ^ see Eq. ( IIO8I 1. 

If the trapping stength increases, one cannot employ pertur- 
bation theory anymore. Nevertheless, one can always calcu- 
late the eigenvalues of H numerically. As it will turn out, the 
scaling 7; ^ still holds in this case and extends even over 
a wider range of ^-values. Figure [32] shows the spectrum of 
7/ for N — 100 and F = 1; the double logarithmic plot (see 
inset) demonstrates that scaling holds for 10 < ^ < 60, where 
the exponent /i is about /i — 1.865. 

Figure[33]compares. for a linear system with N — 100 and 
r = 1, the classical PMit) to the quantum mechanical sur- 
vival probability UMit) |50]. Evidently, Pa/(0 and IlM{t) 
differ strongly: the PM{t) decay established for CTRW is 
practically exponential. llM{t), on the other hand, shows two 
regimes: a power-law decay at intermediate times (panel (a)) 
and an exponential decay (panel (b)) at very long times. 

Turning now to the parameter dependences of 11^/(0' 
Fig.[34lshows the dependence ofIlM{t) on N [50]. Note that 
the scaling regime, where IlM{t) ^ t^^^^ holds, gets larger 
with increasing N . The cross-over to this scaling region from 
the domain of short times occurs around t « A^/2. For larger 
N and in the intermediate time domain, IiM{t) scales nicely 
with . In this case, the power-law approximation [Eq. ( IIO8I 1I 
holds and by rescaling I to l/N one has from Eq. ( I107l i that 
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FIG. 33: (Color online) Temporal decay of Hm (t) (solid black lines) 
and PM{t) (short dashed green lines) for A'^ = 100 and F = 1 in (a) 
double logarithmic scales and in (b) logarithmic scales. Indicated are 
the fits to n M (t) (long dashed lines) in the intermediate (upper red) 
and the long (lower blue) time regime. From 1501 . 



^exp 

I 



2{l/N)f'N-^^-''h 



(124) 



where it was assumed that a ~ N^^ for a Unear system lll09ll . 
Thus, when rescaling / to / /N, the time has to be rescaled by 
the factor N~'^^~^\ Indeed, all curves for which a power- 
law behavior is visible fall on a master curve; see the inset in 
Fig. [31 



3. Line with traps and long-range interactions 

As mentioned above, the interaction range does not need to 
be restricted to nearest neighbor interactions (NNI). When the 
interactions between two nodes go as R^'^, the Hamilto nian of 
the networks without traps has the following structure lll05ll : 



N 



n — 1 

YR-''{\n){n\~\n-R){n\ 



N- 



Y R'''{\n){n\ - \n + R)(n 



(125) 



Note that in the case of a line, the states \k) with k ^ 
{1, . . . , are implicitly excluded from the summation. 

For large exponents v the LRI can be regarded as a small 
perturbation to the NNI, i.e., having Ho{i') ~ Hq + Hy, 




rescaled time t/f^'^ 

FIG. 34: (Color online) Panel (a) shows the A'^-dependence of YVm (t) 
for r = 1; A'^ increases in steps of 10 from 20 (blue line) to 
100 (green line). Panel (b) shows IlMit) versus the rescaled time 
t/iV^-". From (H. 



where Hi, contains only the correction terms to the NNI case 
Hq. This allows one to calculate from the unperturbed states 
the perturbed eigenstates l^*/) up to first order. Taking 
the states to be the eigenstates of the LRI system with- 
out traps, one readily obtains the imaginary parts ji for small 
trapping strength from Eq. (1103b as 7/ = 2r| | , where 

li05D 



(i|*;) = (i|*. 



(0)\ 



E 



p(0) pv 



(0) 



(0)\ 



(126) 



For large v the coupling to the next-next-nearest neighbor is 
by a factor of (3/2)"^ smaller, for = 10 this is about one and 
a half orders of magnitude. Taking, for fixed v, only nearest 
and next-nearest neighbor couplings into account allows one 
to obtain simple analytic expressions. Thus, Eq. ( I126I I yields 



+2-^y^sin(20,)sin(|), (127) 
(6*/ = 7r(iV - e [0, 7r[) which results in |[lol 



(128) 



where 7^ '^^ is the NNI expression, discussed in Sec. |VI B 21 

and 7,'^^ = (8r/7V) cos (6l,/2) sin (261/) sin (6l,/2) is the cor- 
rection due to the LRI. The smallest 7; -values are those for 
which I <C N, which leads to a decrease of the imaginary 
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(a) r=0.00l 





FIG. 36: Examples of two regular hyperbranched fractals of genera- 
tion g — 3 with functionalities f — 3 (left) and / = 4 (right). 



N = 100 five nodes randomly as trap nodes, Xu has shown 
numerically that both the CTQW an d the CTRW survival 
probabilities decay with increasing m ill lOll . thus showing an 
opposite effect to that mentioned above. However, an expla- 
nation was not given. 



FIG. 35: (Color online) //-dependence of the quantum mechanical 
IlM{t) and the classical PM(t) decay behaviors for a chain of A'^ = 
100 sites; here (a) T = 0.001 and (b) T = 1. The inset in (a) shows 
a close-up picture of the region where nAf(t) and PKiit) cross. The 
inset in (b) shows power-law fits to YlM{t) in the intermediate time 
regi me wi th expo nents l//i, where the are taken from Fig. 3(b) of 
Ref. (lOl. From fiol. 



parts 7i because 7^'' < for / <Si N . Here, one can approx- 
imate the imaginary parts by a power-law, i.e., 7; ~ Z^. A 
rough estimate of the scaling exponent /x, assuming v ':$> 1, 
can be readily given IllOSil : 



11: 



ln7;-|-i - In 7; 
ln(Z + l)-lnZ 



(o)+2-V^. 



(129) 



Since ^S^^ is strictly positive for small /, the inclusion of LRI 
leads to a decrease of 7; when compared to the NNI case. In 
turn, this results in a slower decay of 11 a/ (t). 

Figure[35]displays for comparison the quantum mechanical 
IiM{t) and the classical PM{t) behaviors for different ly and 
F; IlM{t) and PM{t) were obtained by numerically diagonal- 
izing the corresponding Hamiltonian H{u) and transfer ma- 
trix T(j/), respectively |105]. Clearly, for both F-values the 
LRI lead to a slower decay of IIm (f), i.e., to a slower trapping 
of the excitation, which is somewhat counterintuitive since the 
opposite effect is observable for classical systems, where the 
decay of Pm {t) becomes faster for decreasing see below. 
By increasing the trapping strength F, the difference between 
the quantum and the classical behaviors become even more 
pronounced, compare Figs. [35]' a) and ISSlb). Generally, for 
nj\/ (<) a change in F leads mainly to in a rescaled time axis, 
since the imaginary parts 7; turn out to be of the same order 
of magnitude when rescaled by F. 

Xu has studied a ring wher e the interaction strength does 
not decay with distance, as in llOSll . but ranges from a given 
node to the m nearest neighbors jl IQl . By choosing for 



C. Fractals 

1. Hyperbranched fractals 

Unlike the regular network discussed in the previous sec- 
tions, hyperbranched fractals allow one to study the survival 
probability in the presence of highly degenerate eigenvalues. 
Volta investigated so-called regular hyperbranched fractals, 
see Fig. [36] for which the eigenvalue spectr a of the connec- 
tivity matrices can be calculated recursively ll 1 ill . As shown 
previously [112], the eigenvalues of the {g + l)st generation 
can be obtained from the eigenvalues of the gth generation by 
solving 



(130) 



where P(A) = A(A — 3) (A — / — 1), / being the functionality 
of the fractal. The thre e roo ts of the polynomial are given by 
the Cardano solutions ll 12ll 

A, = ^ + ^|/(/-l) + 7|i/2cos[(^ + 2^j)/3], (131) 

with j = 1, 2, 3. Thus, each eigenvalue of the gth generation 
gives rise to three new eigenvalues, not all of which are dif- 
ferent from the previous ones. Therefore, there appear highly 
degenerate eigenvalues, e.g., the eigenvalue A = 1 has in the 
gth generation the degeneracy = (/ — 2)(/ + 1)^^^ + 1. 

Placing now a trap node at the center of the fractal, Volta is 
able to calculate the (complex) spectrum of the Hamiltonian 
including the trap 1 1 1 1]. It turns out that only the nondegener- 
ate eigenvalues depend on the trapping strength F. Thus, the 
degenerate eigenvalues can be again calculated based on the 
Cardano solution. In particular, these eigenvalues are real and, 
therefore, do not contribute to the decay of the survival prob- 
abihty. At generation g there are (3^^ - l)/2 + (/ + 1)^* - 3^* 
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real degenerate eigenvalues, such that Eq. ( I107l i yields 

HmW « ^[(3^-l)/2+(/ + l)3-3« 



■ ^ exp(-27,t) 



(132) 



Clearly, in the limit of t cx) this leads to a constant value. 
As Volta also shows (extending the results for the line), in- 
creasing the trapping strength above the value of F = y = 1 
does not lead to a faster 11^/(0 decay. 



D. Random networks 
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position 
of trap 

••t . 



# V \* • • 



1. Disordered system with one trap 

Introducing long-range interactions also allows one to study 
topologically disordered systems. To illustrate this, take a ran- 
dom confi guration of [N — 1) identical nodes and one trap 
node flOT'l, see Fig. [37] All N nodes are placed at random in 
a 3-dimensional box with Cartesian coordinates {x^^^}, with 
i = 1,2, 3. Then the distance between two nodes j and k is 
given by 



3 



-I 1/2 



where the coordinates x\ and x 



ii) 



(133) 



n(0 



"J 

tributed random numbers in the interval [0, N 



are homogeneously dis- 
To relate this 

to the energy transfer dynamics within Rydberg gases one 
considers interactions decaying as RJ'^- In the absence of 
traps the corresponding Hamiltonian Hq has the following 
matrix elements 



R 



{k\Ho\j) 



E 



for k ^ j 



R'l fovk^j. 



'j,k 



(134) 



Now, choose one of the N nodes to be a trap, i.e., for this 
node the full Hamiltonian H has an additional purely imag- 
inary matrix element — iF. Since the configuration of nodes 
is random, one can (without any loss of generality) assume in 
the following that the node labeled 1 is the trap. 

For disordered systems, one calculates averages over TZ dif- 
ferent realizations following again, as for the SWN, 



(135) 



where [ • ] ^ denotes the realization r. In the calculations one 
also assumes the trapping strength [F] ^ to be realization de- 
pendent, because it is required that it be proportional to the 
diagonal element [(l|i?o|l)] in that particular realization, 
namely [F]^ = F [(l|i?o|l)] ^- As it turns out, the depen- 
dence of the decay on the value of F i s quit e weak - different 
F mainly rescale the time axis. In Ref. llOTll only two extreme 



FIG. 37: Random configuration of 99 nodes and a single trap. From 
flOTll. 



cases were considered: (a) F = 10~^, for which a perturba- 
tion theoretical treatment can be justified, and (b) F = 1, such 
that the average trapping strength is of the same order as the 
diagonal elements of Hq at the node of the trap. 

While each realization [n(t)] ^ leads to a specific spectrum 
of the ['yi\^, the relation between {n{t))'jz and the average 
{li)n is not that straightforward. However, for all t, the func- 
tion exp(— 27;i) is convex, therefore, Jensen's inequality ap- 
plies, s ee pa ragraph 12.41 of JtsIi . such that one obtains for a 
given I (im 



(exp(-27ji))^>exp(-2i(7,)fl). 
From this one gets a lower bound for (n(i))^: 



1 ^ 



(136) 



(137) 



In Fig. [38] {n{t))Tz is displayed for the two values F — 
10~^ and F = 1, for different N. In all cases the intermediate 
time decay can be fitted by a power-law 1 10711 

(n(<))K^<-''W^ (138) 

where, different from the regular linear case ifsol llOSll . the ex- 
ponent T] is now iV-dependent. It turns out that approximating 
the exponent by 



(139) 



while keeping F fixed, reproduces the curves well. Note, how- 
ever, that rjo and also /i can still depend on F. 

Now, one can estimate fj, based on the results of rj for 
Ni = 100 and N2 = 1000. Based on ?/o = r]{Ni)/N^ = 
ri{N2)/N^ one has HotIi 



In7/(7V2) - lnry(7Vi) 



(140) 



In N2 - In Ni 

From the numerical values given in Fig. [38] one obtains ap- 
proximately /i w —0.166 and rjQ — 0.0349 for F = 10^^ and 
/i w -0.204 and 770 = 0.0313 for F = 1. 
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FIG. 38: (Color online) Ensemble averages (n(f)) r for different A^: 
(a) for r — 10^® and (b) for F = 1. The scalings in the intermediate 
time regions for A'^ — 100 and A'^ = 1000 are shown as solid black 
lines along with the appropriate scaling law. The arrows are guides 
to th e eye pointing at the bend of (n(t))_R for = 1000. From 

fiol. 



when integrating along the if -axis one has, llll4lll5l] . 

dKW{X,K;t) = \i;{X;t)\^. (142) 



In the case of a discrete system, given, for instance, by N 
discrete positions on a network (enumerated as 0, 1, . . . , — 
1) the functions ip{x) are only defined for integer values of 
X = 0, 1, . . . , - 1, and the form of Eq. (ll41b has to be 
changed from an integral to a sum. Th ere have b een several 
attempts in doing so, see, for instance, llll6llll7ll . However, 
the definition of the discrete WFs might depend on whether 
the length N of the system is even or odd [118, 11^. 

For a one-dimensional system of length N with periodic 
boundary conditions (exemplified by a ring) one has ^"(2;) = 
4>{x zt rN) for all r e N. It follows that each and every 
one of the products ip*{x — y';t)ip{x + y';t) is identical to 
(at least) one of the N forms ijj*{x — y; t)ip{x + y; t), where 
y = 0,l,...,N-l. 

Now the WF has the form of a Fourier transform; a unique 
transformation of these N products requires N different k- 
values. These fc-values may evidently be chosen as fc = 
2ttk/N, again having k = 0, 1, . . . , — 1. One is thus led to 
propose for integer x and y the following discrete WF Iil20il 



y=o 



From Fig. [38] as well as from Eq. ( |139t one sees that the 
exponent 77 (iV) decreases with increasing N. Certainly, if N 
becomes very large it becomes quite improbable (in the en- 
semble average) for an exciton to encounter the single trap. 
Therefore, the decay of (n(<))/f can only be observed at very 
long times. 



VII. RELATIONS BETWEEN CTRW/CTQW AND OTHER 
APPROACHES 

A. Phase space approaches 

In order to obtain a unifying framework for describing both 
quantum mechanical and classical transport one introd uces a 
quantum analog for the classical dynamics, see, e.g. ll 1 1 3[l . 
One particular approach is attributed to Wigner and uses the 
so-called Wigner-function (WF) in the (quantum mechanical) 
2d phase space. If the phase space is spannedj^y the c ontinu- 



ous variables X and K, the WF is given by 



id by the c c 



W{X,K;t) J dY , 



iKY 



{X-Y/2\p{t)\X + Y/2), 
(141) 

where p{t) is the density operator and thus for a pure state, 
p{t) = \'tp{t)){ij{t)\. Here, ij{X;t) = {X\i:{t)) is the wave 
function of the particle. The WF is a quasi-probability (in the 
sense that it can become negative). Integrating W{X, K; t) 
along lines in phase space gives marginal distributions, e.g.. 



1. WF for a ring 

Since the eigenstates of the ring are Bloch states, the WF 
for a CTQW on a ring of N nodes reads il2Ql 



N-l 



{x,K;t) = ■^Y.exp[-t27T{2n + K){x-j)/N] 



71=0 



X exp { — i2i[cos(27r(K + n)/N) 
-cos(27m/7V)]}, (144) 

where it has been used that 6 — 2TTn/N and k = 2tik,/N . 
Furthermore, one may note from Eq. ( 11431 ) that 

Y,Wj{x,n-t) = \^{x;t)\^ 

K 

^ E E e-'^^^y^^'r {x - y; mx + y; t) 



(145) 



In general, the WFs have a very complex structure. Figure 
|39] shows a contour plot of the WF of a CTQW on a cycle of 
length N = 101 at different times. Note that at f = the WF 
is localized on the strip at the initial point j. At < = 1, the 
WF is still mostly localized about j. As time increases other 
sites get populated. On short time scales, the WF develops a 
very regular structure in phase space, with "wavefronts" orig- 
inating from the initial point j. Additionally, one also notes 
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x=50 ^ initial node 




FIG. 39: (Color online) WFs of a CTQW on a cycle of length iV = 
101 at times t = 1,10,20,40 [(a)-(d)] as well as t = 100,500 
[(e), (f)]. The initial node is at j = 50. Red regions denote positive 
values of Wj {x, k; t\ blue regions negative values and white regions 
values close to 0. From |[12Q1 . 
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FIG. 40: (Color online) Same as Fig[39l for N = 100 and j = 50. 
From 112(1 . 



"wavefronts" starting from the region opposite to the initial 
point, which are much weaker in amplitude. As time pro- 
gresses, these two types of waves start interfering with each 
other. 

FiguregOlshows the WF of a CTQW for iV = 100. At short 
times, the structure of the WF is quite similar to the situation 
for iV = 101. Nonetheless, there are differences at larger 
times, visible by comparing Figsl39|e) and [39l f) to Figsl^Ole) 
andgaf). 

However, although at long times the interference effects are 
quite intricate, typical patterns are still visible. Kvt — 500, 
one finds less regularities in phase space for N — 101 than 
for N = 100, reflecting the higher symmetry of CTQW for 
even N . One notes, moreover, that the phase space patterns 
give a much richer picture of the underlying dynamics than 
the transition probabilities (a;; f)p alone. 

Furthermore, one should remark that in an infinite system 
the WFs have a much simpler structure, because there one 
does not face the problem of distinct wave fronts running into 
opposite directions and interfering with each other because of 
the closure of the ring. 

In general, the marginal distribution Wj (x, k; t) is ob- 



tained from Eq.(fT44ll as lfT20ll 



W,{x, K;t)^j^J2 exp[z2^n(x - 3)/N] 

K n—Q 

X exp[~i2t 005(27:12/ N)] 
1 

X — ^ exp[i27r(K + n){x- j)/N] 



X exp[i2t cos(27r(K + n) /N)] 



(146) 



Since the system is periodic, the arguments (k + n) in the 
exponentials can be written as (k + n) = N — v, where 
= 0,1,..., TV - 1. This yields 



1 ^-1 

= -y 



gi27rn(2;-j)/JVg-i2t cos(27rn/JV) 



(147) 



which is exactly what also follows, see Eq.( il45l l, from calcu- 
lating_k/;(a;; i)p directly from the Bloch ansatz, see Eq.lfTTb 
and 117211 . 

For the marginal distribution Wj{x, k; t) one also gets 
fromEq. (ll44b 



Af-l 



^i27r{2n+K)j/N 



Tl = 



X exp { - i2i[cos(27r(K + n)/N) - cos(27m/A^)]} 
1/iV for odd and aU K 

2 /N for A^ even and n even (148) 
for A^ even and k odd, 



= < 



all of these expressions are independent of t. Eq. (1148b can 
be confirmed directly by taking the Fourier transform of 
\ip{x; The whole phase space volume is normalized to 
unity, as can be seen by summing Eq. ( 1148b over all k, with 
K = 0, 1,...,A^-1. 

In the long-time average and for odd A^ (superscript °) most 
points in the quantum mechanical phase space have a weight 
of 1 /N^, namely one has has Ill03ll 



1/A^^ for K ^ and any x 
W°{x, k) = {l/N for K = and x = j 
else. 



(149) 



For even A^ (superscript ^), the limiting WF reads 

K 7^ 0, K even and any x 
l/N K = and a; = j, j + N/2 (150) 
else. 

The long time averages of the WFs for even A^ are somewhat 
peculiar, since values different from zero appear only for even 
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K, whereas the WFs themselves have values different from 
zero at arbitrary times for all k. These stripes in the long-time 
average are due to the periodicity of the ring. For even N 
one finds constructive interference patterns in the transition 
probabilities, since the number of steps in both directions is 
the same, see also Ref. ll72ll . For a finite line with even N, 
there are no stripes in the long time average. 



2. Rings with energetic disorder 

By introducing (static) disorder i nto the system, as in 
Sec. IV B[ the Bloch property is lost flo5l. In order to have 
a global picture of the effect of the disorder on the dynam- 
ics, one considers ensemble averages of the WFs. For this one 
calculates the WF for different realizations of H and averages 
over all realizations, i.e., for R realizations: 

1 ^ 

{Wjix,n;t))n ^ - ^ (151) 

r—l 

where [Wj {x, K;t)\^ is the WF of the j'th realization of H. 

It turns out that the particular type of disorder, i.e., diago- 
nal or diagonal and off-diagonal, does not change the picture 
significantly. For diagonal and for off-diagonal disorder Fig- 
ure |4T] shows snapshots of {Wj{x,K;t))fi at different times 
for = 101 and for different values of A. 

The first column shows {Wj {x , k; t)) for A = 1/40 at 
times i = 1,10,20,40,100, and 500 [Fig. |4nia)-(lf)]. For 
this quite weak disorder, the patterns in phase space are sim- 
ilar to the unperturbed case, where "waves" in phase space 
emanate from the initial site x = j = 50 and start to in- 
terfere af ter h aving reached the opposite site of the ring (see 
Fig. 3 of II120II ). However, at longer times differences become 
visible. Fig ITlTlf). The pattern for the unperturbed case is 
quite irregular but with alternating positive and negative re- 
gions of the WF of approximately the same magnitude. For x 
close to the initial site j = 50, the disorder causes a decrease 
of {Wj{x, k;<))_r for K-values in the middle of the interval 
[0, — 1] when compared to what is found for values of k 
close to or to — 1. 

Increasing the disorder parameter A, the patterns change 
profoundly. The wave structure gets suppressed and for all k 
a localized region forms about the initial site j = 50, already 
for small disorder (A = 1/10) and short times {t = 20), see 
Fig.|4ll2c). 

For even larger values of A, for all k the formation of a lo- 
calized region about j = 50 becomes even more pronounced. 
Already for A = 1/4 this localized region forms for times as 
short as t = 10, see Fig.|4n3b). At A = 1/2, {Wj{x, k; t))ii 
stays localized at all times [Fig. |4Tj4a)-(4f)]. Also here, the 
values of the WF at about k w N/2 are rather low, whereas 
the values of the WF for k close to the interval borders and 

— 1 remain rather large for x « j, as indicated by the thin 
black region (see also the arrows), e.g., in Fig. 141140. One 
further notes that at high disorder the localized averaged WF 
is always positive, i.e., all fluctuations, present for small dis- 
order, have vanished. One recalls that the WF is normalized to 
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FIG. 41: (Color online) Ensemble average of WFs of the quantum 
dynamics on a ring of length A'^ — 101 with diagonal and off- 
diagonal disorder for A = 1/40, 1/10, 1/4, and 1/2 [columns (1)- 
(4)], each at times i = 1, 10, 20, 40, 100, and 500 [rows (a)-(f)]. The 
initial node is always j — 50 and the average is over R — 1000 re- 
alizations. Red regions denote positive values of the averaged WFs, 
blue regions negative values and white regions values close to 0. The 
colormaps are always chosen to be the same for each row but might 
differ in different rows. The maximal values of {Wj{x, k; f))_R are 
denoted by small black regions; these are highlighted and exempli- 
fied by the arrows in panel (4f). From 11031. 



unity when integrated over the whole phase space. Having an 
even number of nodes in the graph does not alter the picture 
significantly. 

The ensemble average of the long-time averaged WF fol- 
lows as ifTol 



{Wj{x, k))r = (\im ^ [ dtWj{x,K-t)] 



lim — 

T~>oo T 



dt {W,{x,K;t))i 



(152) 



Now, the disorder changes also the limiting WF quite drasti- 
cally. Starting from high disorder of A = 1/2, one expects 
from Figs. l4TT 4a)-(4f) that the long time average of the av- 
eraged WF will look basically the same. Figure |42] shows 
the limiting averaged WF for A^ — 101 and DOD accord- 
ing to Eq. (|152| i. Indeed, for large A, the limiting averaged 
WF is comparable to the corresponding averaged WF, com- 
pare Figs. I42l 4). Close to the initial node a; = j = 50 and 
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FIG. 42: (Color online) DOD: Limiting averaged WF, {Wj{x, k))r, 
for N = 101 and A = 1/40, 1/10, 1/4, andJ/2 [panels (1) to (4), 
respectively], according to Ea. U52t . From 110311 . 



see Eq. dlSOI l. By switching on the disorder, these peculiar- 
ities vanish. Figure |43] shows the Hmiting averaged WF for 
DOD and N = 100. Ahhough for small A there are some re- 
mainders of stripes left, these disappear completely for higher 
values of A, compare Figs. l43l' l)-(4). Note further that the 
second peak of Wj{x,K) at x = j + N/2, see Eg. ( 11501) . 
transforms for increasing disorder to an oscillatory line in the 
K-direction at a; = j + N/2. 



3. Long-range interaction cycles 

X u and Liu have considered the WF, defined in Eq. (11431 ) 
and II120II . for rings, where the m next nearest neighbors of a 
given node are connected to this node by bonds lll2lll . see also 
Sec. IV Al The eigenstates of such a translationally invariant 
system are still Bloch states. However, the eigenvalues read 
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FIG. 43: (Color online) Same as Fig.|42]but for = 100. From 
{M\. 



there along the K-direction, {Wj{x, k))r has large (positive) 
values for k about and — 1 which decrease by going to- 
ward K = N/2. Here, the decrease depends on the particular 
type of disorder involved. 

Also for small A there are significant differences to the 
case without disorder From Fig. \42\ l) one sees that the 
disorder "smears out" the localized value Wj{j,0) — 1/N 
in the absence of disorder, see Eq.( |149t . Specifically, the 
{x = j,K — 0)-value decreases, while the neighboring ones 
increase. For A = 1/10, the onset of localization about 
X = j — 50 can already be seen, e.g., Fig. l42l' 2). and becomes 
more and more pronounced as A increases. Fig. l42l" 3). Fur- 
thermore, all other values of {Wj {x,k)) r for x ^ j decrease 
with increasing disorder, as can be seen by the decreasing size 
of the light pink region, which corresponds to values close to 
1/iV^. In fact, the values of the limiting averaged WF for x 
distant from a; = j = 50 drop to zero, shown as white regions 
inFig.l42l 

For even one finds that without disorder the limiting WF 
shows a peculiar "striped" distribution, caused by the PRCs, 



^„ = 2m - 2 ^ cos 2nji:/N. 



(153) 



As in Eq. (I144l i. the WF follows as 



N-l 



W,{x,.:t) = 1 ^e-*(^"-^~- 



iV2 

xe 



n=0 

i2ir(2n-K)(x-j)/N 



(154) 



Thus, inserting the results for the eigenvalues one obtains the 
WF for this type of long-range interacting cycles. 

Now, while for ra — 1 one recovers the results of the pre- 
vious section, increasing m leads to more complex structures 
of the WF. The fronts emanating from the initial value x — j 
have, in the case of to = 1, only a single maximum value 
at about k — N/2, see Figs. |39]andl40l For m > 1 there 
appear to fro nt ma xima at about k ^ {21 — l)N/2rn with 
I = 1, . . . , m With increasing number m the structure 

of the WFs become more complex, but the resulting patterns 
can stil l be attributed to the interfering fronts emanating from 

x = j iili. 



4. Disordered networks 

Xu and Liu have also studied the WF for disordered sys- 
tems, namely for the Watts -Strogatz (WS) small-world model 
II121II . The WS model starts from an ordered ring with only 
nearest neighbor interactions. Now, every bond is rewired 
with probability p, implying that for p — one retains the 
ordered ring and for p = 1 one has a completely disordered 
network of N nodes with N bonds. The WF for such systems 
shows essentially the same features as for rings with energetic 
disorder, see above. In the ensemble average there is stronger 
and stronger localization with increasing p. Moreover, here 
also the distinction between even and odd N remains, namely 
even N lead to a stripe-patterned WF caused by co nstruc tive 
interference while odd N do not show this behavior lll2lll. 
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B. Quantum master equations 



1. Decoherence on rings 



CTQW and CTRW are the two extreme cases of possi- 
ble transport processes, i.e., a purely coherent transport for 
CTQW and a purely incoherent transport for CTRW. Now, 
if a quantum mechanical system cannot be regarded as being 
isolated from its environment, this environment will influence 
the dynamics of the system. Depending on the strength of the 
coupling (where also the temperature is assumed to be con- 
trolled by the environment), the dynamics may change from 
quantal to classical. Since the total system (i.e., the system 
S and the environment/reservoir R) resides in a Hilbert space 
of huge dimension (this due to the many degrees of freedom 
(DOF) of the environment), one often considers the system 
only, where the environmental DOF have been traced out |69]. 

Let W{t) be the density operator of the total system, then 
p{t) = ti-fj{W{t)} is called the reduced density operator 
Then the dynamics of W{t) is still governed by Schrodinger's 
equation, which can be recast into the Liouville-von Neumann 
(LvN) equation 

Wit) = -t[Ht,t,W{t)], (155) 

where iiftot is the Hamiltonian of the total system. Jftot usu- 
ally comprises three parts: the Hamiltonians Hs for the sys- 
tem, Hii for the environemnt, and Hs for the coupling be- 
tween system and environment: 



Htot — Hs + Hji + HsR- 



(156) 



Inserting this into the LvN equation for W{t) and integrat- 
ing out the environmental DOF one arrives at an equation for 
the reduced density operator (for details see ll69ll ): 



(157) 



Here, 2?[p(t)] is called the dissipator. Note that in order to 
arrive at Eq. ( 11571 ) several assumptions, such as bilinear cou- 
plings between system and environment or the Markov ap- 
proximation, have been used . 

Under certain conditons, such as a weak coupling between 
the system and the environment, the dissipator can also be ex- 
pressed by operators acting on the Hilbert space of the system 
S. In particular, Eq. (11571) can be brought into the so-called 
Lindblad form ifiil 



Pit) = -i[Hs.p{t)]~Y.{^^^P^^)^ 
-L]L,p{t) - p{t)L]L, 



(158) 



where the operators Lj, the Lindblad operators, mimick the 
influence of the environemnt on the dynamics. Equation ( |158t 
is also called Lindblad quantum master equation (LQME). 

Obviously, if the couphng to the environment vanishes, one 
arrives again at the LvN equation for the density operator p{t) 
of the system S. Since the resulting LvN equation is equiv- 
alent to the Schrodinger equation of a closed system with 
Hamiltonian H, it is also equivalent to the CTQW formula- 
tion of the previous sections. Now, when increasing the cou- 
pling to the environment, Eq. ( |158t allows to study the onset 
of the quantum to classical cross-over. 



Take now a ring of N nodes and N Lindblad operators each 
of which projects onto a single node, i.e., Lj = ^/Aj | j) {j \ for 
i = 1, . . . ,N. The parameters Xj specify the strength of the 
coupling to the environment. If Xj — X for all j, Eq. ( 11581 ) 
simplifies: 



Pit) - -i[Hs,p{t)] 

-2XY,{p{t)~p,,{t))\j){j\, 



(159) 



where pkj{t) = {k\p{t)\j). The rate A can be estimated from 
the spectral density J{uj) d escri bing the environment within 
the Caldeira-Leggett model lll22lll23[l at a given temperature 
T. Taking J{u!) = 27raaj exp(— w/wc) and using the Markov 
approximation one arrives at A = TTaksT ll69ll . One has to 
bear in mind that Eq. ( |159t is an approximation with a lim- 
ited range of validity: For a very large coupling strength A, 
Eq. ( 11701 ) leads to the quantum Zeno limit rather than to a clas- 
sical master/rate equation. In matrix form Eq. ( 11591 ) is equiv- 
alent to the s o-called G urvitz model, see for instance Iil24il . 
considered in il25l4l27ll 



Pkjit) - -i{k\[Hs,pit)]\j) 
-2Xil-6kj)pkjit). 



(160) 



Solenov and Fedichkin showed analytically that the prob- 
ability distribution to be at node k, where for a ring one as- 
sumes without loss of generality that the initial node is j = 0, 
foll ows from t he diagonal elements of the reduced density ma- 
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2TTi , , , ■ 

--^(n + TOjfe 



(161) 



Clearly, Eq. (1161b reduces to the results of Sec. IIV Al in the 
limit of A ^ 0. For A 7^ 0, due to the exponentially decaying 
factors, the probability distribution will eventually decay to 
the classical equilibrium values, in this case the equipartition 
value 1 IN. 

The so-called mixing time t,„ is used to compare the 
spreading of CTQW to that of CTRW. It is implicitly defined 
by looking at the deviation from the classical uniform distri- 
bution XjN [125il . One gets thus 



Af-l 

E 

fc=0 



I ^ 1 

Pk,k(tm) - 



<£, 



(162) 



where e is a small dimensionless constant representing the de- 
gree of mixing. By using Eq. (1161b . Solenov and Fedichkin 
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showed that the mixing time is bound from above by 



2. Dimer with traps 



minim < 



N 



2X{N -2) 



In 



iV+ 1 



(163) 



The same authors have generahzed this approach to so- 
called hypercycles. A hypercycle of dimension d and size N 
is build by considering N copies of the d—1 dimension al hy- 
percycle and connecting the corresponding nodes, see lll25ll 
for details. Thereby, the total number of nodes is N'^. The 
corrsponding upper bound for the mixing time tm now fol- 
lows as 



'"''^ - 2A iV - 1 



In 



d{N + + eN'^) 



(164) 



Within the phenonmenological approach presented above, 
introducing traps into the system leads to non-Hermitian 
Hamiltonians. The LvN equation changes, therefore, to 



W{t) = -z[Ho,W{t)] - {T,W{t)}, 



(169) 



where {•, •} is the anti-commutator and Hq the Hamiltonian 
without traps. Writing the L indblad operators in the form 
VXLj, the LQME reads CH 

p{t) = -i[Ho,pit)]-{T,pit)} 

-2\Y,(p{t)-p,,{t))L,. (170) 



Salimi and Radgohar have considered the long-range case 
of the single cycle, where a given n ode i nteracts with its 2m 
nearest neighbors (m to each side) f 128*], see also Sec. IV Al 
The probability distribution to be at node k is then obtained 
for even m as 



N 

X exp 
X exp 
X exp 
and for odd m as 
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(165) 
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While the upper bounds of the mixing time turn out to be in- 
dependent of TO, they do differ for even and for odd to, one 
has namely 



< 



N N 



and 



< 



N 



N 



2 e 



for even to 



for odd TO. 



(167) 



(168) 



Therefore, mixing turns out to be faster for even to than for 
odd TO. 



Now, consider a dimer which is coupled to an external bath. 
This situation allows to solve Eq. ( 11701 ) analytically and to 
compare the LQME results to the numerically exact Path In- 
tegral Monte Cai-lo (PIMC) calc ulations. For details on PIMC 
techniques see, e.g., 1 130. Il3l[l and references therein. The 
Hamiltonian of the dimer without any coupling to the sur- 
roundings can be expressed through the Pauli matrices cTz and 



H ^ El 



VcTx - i-^i^ - CTz) 



(171) 



where E is the on-site energy, which is choosen to be the same 
for both nodes, and V is the coupling betw een th e two nodes. 
It is easily verified that the eigenvalues are lll29ll 



E±=E± Ve^"'' = E± y/V^ - r2/4 - iT/2, (172) 

where (/> — arcsin(r/2F). For F — (0 — >• 0) this yields 
the correct eigenvalues E ± V of Hq. Note that for F < 
2V the negative imaginary part of E± is identical for both 
eigenvalues, i.e., 7+^7-= F /2. T he bi-orthonormalized 
eigenstates of H are of the form lll29ll 



1 



^^^^ ^ \±e^'*/' 



and 



(173) 



(174) 



where the phases (f) depend on F, such that in the limit F — > 
one recovers the eigenstates of Hq. 

When the coupling to the environment vanishes (A — > 0), 
one obtains the survival probability directly from the eigen- 
states and eigenvalues of H. For F < 2V one has [rl29il 



^, , cos^ ((h + tV cos ( 
n(t)=e ^* ^— — — : 



(forA = 0). (175) 



Note that for values F > 2V the dimer is overdamped. 
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When considering the dimer without traps (F = 0) but cou- 
pled to the environment, Eq. jllOi simpHfies and, from the 
solution for p, one obtains the transition probabilities [ 129.1 



(a) 



-nit) 



A sin {WW^ - A2) 



(forr = 0). 



(176) 

Moreover, tt^I (t) = 1 — irf'l (t). From Eq. ( I176I I one recovers 
for A — ^ the simple oscillatory behavior of the transition 
probabilities [namely, \im\^Q Tr^^l{t) = cos^(yt)]. When 
the coupling to the surroundings does not vanish, for A > 0, 
the transition probabilities still show oscillations, which are 
superimposed on an exponential decay which tends for long 
times to the classical equipartition value of 1 /2. 

In order to combine the results of Eq. ( I175l l - taking val- 
ues r < 2V - and Eq. ( 11761 ), one expands in both equations 
all terms except the exponentials to first order in F and A, re- 
spectively. Note that for Eq. ( 11751 ) one obtains a product of 
exp(— Ft) and cos^ Vt, which is the simple oscillatory behav- 
ior of the dimer in the absence of the trap. Now, the coupling 
to the environment affects all transitions but still conserves 
probabilities. Therefore, one replaces the term cos^ Vt by the 
expansion of Eq. ( 11761 ), such that one obtains il29|[ 



-rt (0) 
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■ sin 2Vt 



(177) 



Figure |44] compares the survival probabilities of a dissi- 
pative dimer obtained from the approximative LQME ap- 
proach to the numerically exact PIMC calculations for a bath 
with ohmic spectral density with exponential cutoff, J{lo) = 
2710^6^"/""^. The initial condition is 7ri_i(0) — 1, i.e., at 
t = the system is localized in the non-trap node 1 of the 
dimer Here, the on-site energies E and the coupling elements 
V have been taken to be equal, E — V — 1, one sets ujc ~ 5V 
and the temperature is fixed to fc^T = T^. 

For small trapping strength (F = 0.1) and vanishing cou- 
pling to the environment (a = 0), Fig. |44ta), the PIMC cal- 
culations coincide with the result of Eq. (11751 ). A moderate 
increase of the coupling {a = 1/10), Fig. |44|b), still leads 
for Eqs. (11701 ) (solid lines) and ( 11771 ) (dashed lines) to re- 
sults which are in excellent agreement with the findings of 
the PIMC calculations (symbols). When increasing the cou- 
phng further to a = 1/4, Fig. |44jc), however, the approxi- 
mate solution, Eq. ( |177t , begins to deviate from the LQME 
and the PIMC calculations, which are still in very good agree- 
ment lfT29t[ . 

As the numerical effort of real-time PIMC simulations 
grows exponentially with time, they can cover only short-to- 
intermediate time scales. However, the agreement between 
the LQME and the PIMC calculations in the weak coupling 
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FIG. 44: (Color online) PIMC results (circles) for a dimer with 
F = 0.1 and different system-bath couplings a — X/n: (z) a — 0, 
(b) Of = 1/10, and (c) q = 1/4. The solid lines represent the numer- 
ical solution of the LQME equation, the dashed lines show the corre- 
sponding analytical results obtained from Eq. J175t for a = X = 
and from Eq. ( 1177b for a = 1/10 and a = 1/4. The dotted blue line 
shows the long time limit Il{t) ~ exp{—Tt). Panel (d) shows the 
corresponding long-time behavior of the numerical LQME solution 
for the three different values o f a a nd additionally the behavior for 
large couplings a = 10. From Il29ll . 



regime permits to compensate for this shortcoming; it allows 
namely to use for longer times the LQME procedure, see 

Fig.nad). 



VIII. OUTLOOK 

The application of CTQW to transport for large classes of 
(physical, chemical, and biological) phenomena involving dif- 
ferent types of networks has turned out to be very successful 
in recent years. However, only little is known about the de- 
tailed relations between topology, interaction ranges, dimen- 
sions and the transport efficiencies. While CTRW fall into dif- 
ferent universality classes depending on the (fractal) dimen- 
sions and the interaction ranges, there are only few indica- 
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tions on the analogous emergence of "quantum universality 
classes". Take as an example the dynamics on a chain with 
long-range interactions decaying as \k — j\ see Sec. IV Al 
in this case the mean square displacement of CTQW increases 
as also for 7 = 2. This is not the case for CTRW, where 
only systems with 7 > 3 lead to the same power-law behavior 
of the mean square displacement. 

Therefore, a thorough investigation of the influence of dif- 
ferent topological aspects on the dynamics is clearly neces- 
sary. Moreover, once a systematic classification of CTQW 
gets established, one needs to examine the transition from the 
quantum to the classical picture, in order to understand the 
changes in the universality patterns. For this, the environ- 
ment of each quantum system has to be properly taken into 
account. In a phenomenological ansatz one can incorporate 
this either by using the quantum master equation approach, 
sketched in Sec. lVIIBlo r by using so-called generalized mas- 
ter equations (GME) f^. Here, the Markovian master equa- 
tion (as used in CTRW) is extended by a memory -kernel (MK) 
to account for non-Markovian effects. In the limit of a MK of 
delta-type one recovers the CTRW, while in the limit of a con- 
stant MK one arrives for CTQW at a wave equation similar 
to the Schrodinger equation. Both approaches allow to study 



the influence of the environment on the dynamics. Clearly, 
the validity of both methods has to be checked. One option to 
do so is to compare the results - for not too large systems - to 
those of numerically exact PIMC simulations. 

Such a comparison has several advantages. First, it sets the 
range of validity of the phenomenological approaches. Sec- 
ond, once this range is established, it allows to extend the 
short- time PIMC results to (in principle) arbitrary long times. 
Thus, joining different methods together may lead to a deeper 
understanding of the complex transport processes found in na- 
ture. 
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